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Abstract. We study couplings q' of two equivariant random measures A* and fi' on a Riemannian 
manifold (M,d,m). Given a cost function we ask for minimizers of the mean transportation cost per 
volume. In case the minimal/optimal cost is finite and A" <C m we prove that there is a unique equi- 
variant coupling minimizing the mean transportation cost per volume. Moreover, the optimal coupling 
is induced by a transportation map, q' = (id, T)*A*. We show that the optimal transportation map 
can be approximated by solutions to classical optimal transportation problems on bounded regions. In 
case of L p — cost the optimal transportation cost per volume defines a metric on the space of equivariant 
random measure with unit intensity. 



1. Introduction and Statement of Main Results 

Let {M, d, m) be a connected smooth non-compact Riemannian manifold with Riemannian distance d, and 
Riemannian volume m. Assume that there is a group G of isometries of M acting properly discontinuously, 
cocompactly and freely on M. A random measure A* on M is a measure valued random variable modeled 
on some probability space (0,21, P). We assume that the probability space admits a measurable flow 
(O g ) ge Q which we interpret as the action of G on the support of A w . A random measure A* is called 
equivariant if 

\ e ^( g .) = A w (.) for all u e 0, g e G. 
We will assume that P is stationary, that is P is invariant under the flow 9. In particular, this implies 

that \*{B) = X'(gB) for any g £ G and Borel set B. All random measures will be defined on the same 
probability space. 

We want to extend the theory of optimal transportation to the case of equivariant random measure A*, ji* 
on M. Due to the almost sure infinite mass of and fi u the usual notion of optimality, namely being 
a minimizer of the total transportation cost, is not meaningful. Therefore, we restrict our investigation 
to the case of equivariant random measures. For, equivariance allows to transform local quantities into 
global quantities. To be more precise, if an equivariant coupling can be locally improved it can also be 
globally improved. 

Hence, given two equivariant random measures (A*,//*) of equal intensity on M, we are interested in 
couplings q* of A* and fi' , i.e. measure valued random variables ui H ► q u such that for any u € f2 the 
measure on M x M is a coupling of A" and fj, u . We look for minimizers of the mean transportation 
cost 



€(q') := sup -I-E 

BeAdm(Af) m(tf) 



c{x,y) q'{dx,dy) 

MxB 



where Adm(M) is the set of all bou nded Borel sets that can be written as the union of "translates" of 



fundamental regions (see section 2.8 ). For example for M = K d , G = Z d acting by translation, a typical set 
would be a finite union of unit cubes. We always consider cost functions of the form c(x,y) = $(d(x,y)) 
for some continuous strictly increasing function ■& : M + — > M.+ with #(0) = and linv^oo #(r) = oo. 
Additionally, we assume that the classical Monge problem between two compactly supported probability 
measures A and fi with A <ti m has a unique solution. 

A coupling q* of A* and fi* is called optimal if it is equivariant and minimizes the mean transportation 
cost among all equivariant couplings. The set of all equivariant couplings between A* and \i* will be 
denoted by II e (A*,/Lt*). We will show that there always is at least one optimal coupling as soon as the 
optimal mean transportation cost is finite. A natural question is in which cases can we say more about 
the optimal coupling? When is it unique? Is it possible to construct it? Can we say something about its 
geometry? The first main result states 
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Theorem 1.1. Let (A*, jj,') be two equivariant random measures on M. If the optimal mean transportation 
cost is finite 

c e ,oo = inf < oo 

and A" is absolutely continuous to the volume measure m for almost all oj, then there is a unique optimal 
coupling q' between A* and fx*. It can be represented as q u = (id,T ul )^X u for some measurable map 
: supp(X UJ ) — , supp(p u ) measurably only dependent on the a— algebra generated by (X*,p*). 

In particular, considering A* = m being the Riemannian volume measure and fi* a point process on 
M the optimal transportation map T u defines a fair factor allocation for /z* . The inverse map of T u 
assigns to each point ("center") £ of p? a set ("cell") of volume If the point process is simple, all 

the cells will have mass one. In the case of M = M. d and quadratic cost c(x,y) = \x — y\ 2 all cells will 
be convex polytopes of volume one, they constitute a Laguerre tessellation (see |LZ08| ). In the case of 
linear cost c(x,y) — \x — y\ all cells will be starlike with respect to their center, the allocation becomes 
a Johnson-Mehl diagram (see |AHA92j ). In the light of these results one might interpret the optimal 
coupling as a generalized tessellation. If p' is even invariant under the action of R d the optimal cost 
between the Lebesgue measure L and p* is given by 

Cei00 = E[0(|r(O)|)], 

recovering a quality factor studied by Peres et alii in the context of allocations (e.g. |CPPRb| ). 

Moreover, we prove that the optimal coupling Q°° , if it is unique, can be obtained as the limit of classical 
optimal couplings of A* and p* restricted to bounded sets. For the construction we need to impose an 
additional growth assumption on the group G. The assumptions on the group action imply that G is 
finitely generated. Let S be a generating set und consider the Cayley graph of G with respect to S, 
A(G, S). Let A r denote the closed 2 r neighbourhood of the identity of A(G, S). We will assume that G 
satisfies some strong kind of amenability or otherwise said a certain growth condition, namely 

^ — m — = °' 

r-voo |A r | 

for all g £ G, where | • | denotes the cardinality and A the symmetric difference. Let B be a fundamental 
region and B r = A r B . Let Qb t be the unique optimal semicoupling between A* and that is the 

unique optimal coupling between p ■ A* and 1b,. A** for some optimal choice of density p. Put 



Theorem 1.2. Let (A*,/i*) be two equivariant random measures on M, such that the optimal mean 
transportation cost are finite, c ejOC < oo. Assume, that X u is absolutely continuous to the volume measure 
m for almost all u. Then, for every g £ G 

Q r g -t Q°° vaguely 

inM(M x M x fi). 

For the proof of this theorem the assumption of absolute continuity is only needed to ensure uniqueness of 
Q g B r and Q°° . If we do not have absolute continuity but uniqueness of Q g B r and Q°° the same theorem 
with the same proof holds. 

In the case of absolute continuity we can even say a bit more and get rid of the mixing. The unique 
optimal coupling is given by a map, that is 

Q°° = (id,T)*A*. 
Moreover, the optimal semicoupling Q g B r is given by 

QgB r = (id,Tg,r)*{pg,ry), 

for some measurable map T g ^ r and some density p g , r - Then, we have 
Theorem 1.3. For every g E G 



Tg t r — > T locally in A* (£> P measure 
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Analogous results will be obtained in the more general case of optimal semicouplings between A* and /x* 



where A* has intensity one and /x* has intensity f3 £ (0,oo) (see Theorem 4.11 Theorem 5.3 Proposition 
5.7 and section [6]). In the case (3 < 1, A* is allowed to not transport all of its mass. There will be some 



areas from which nothing is transported and the fi* mass can choose its favorite A* mass. In the case 
P > 1 the situation is the opposite. There is too much \i° mass. Hence, A* can choose its favorite /x* 
mass and some part of the \x* mass will not be satisfied, that is they will not get enough or even any of 
the A* mass. 

As a special case of our result, we recover the results by Huesmann and Sturm in |HS1 . They studied 
couplings between the Lebesgue measure and an equivariant point process of intensity j3 £ (0, 1]. They 
showed that there is a unique optimal semicoupling and also proved an approximation result by solutions 
to transport problems on bounded regions. 

Furthermore, in |HS10j necessary and sufficient conditions have been derived implying the finiteness of the 
mean asymptotic transportation cost in the case of transporting the Lebesgue measure to a Poisson point 
process. By applying the same techniques similar estimates can be achieved for the case of a compound 
Poisson process with iid weights (X^gN, i.e. /i w = X^eN^^i- ^ n ^ ne case * na ^ -^i ^ s exponentially 
distributed it is possible to use the algorithm by Marko and Timar [MTllj to construct an equivariant 
coupling between the Lebesgue measure and /i* with optimal tail, i.e. with finite mean transportation 
cost for the cost function c(x, y) = exp(ft • \x — y\ d ) for some positive k and dimension d > 3. 

In the case of c(x, y) = d p (x, y) for p £ [1, oo) we write the optimal mean transportation cost between A* 
and/x* as W£(A*,/z*), i.e. 

W£(AV) = inf 

Let 

V p = {equivariant random measures p.* on M with unit intensity s.t. W p (m, /j,') < oo} 
Then W p defines a metric on V p which implies the vague convergence of the Campbell measures (see 



Propositions 8.1 and 8.3). Moreover, if we take two sequences of equivariant random measures (A*)„ s n 
and (/i* )„gn such that their Campbell measures converge vaguely to some equivariant random measures 
X',fJ,', i.e. 

A'P -> A'P, ny -> /u'P, 

and the optimal mean transportation cost converge W p (A*,/x*) — > W p (A*,/x*), then the optimal semi- 
couplings between A* and \i* n converge to the optimal semicoupling between A* and \i* (see Proposition 



8.5) 



It is clear that if the restriction of A w -C m is lifted there will not be a unique semicoupling in general. 
However, in the case of two independent Poisson processes on M rf we conjecture that there is a unique 
optimal semicoupling given that the mean transportation cost is finite. In particular, this would imply 
that the optimal coupling between two Poisson processes is a matching. 

Matchings of two independent Poisson processes have been intensely studied in [Hol09, HPPS09J. How- 
ever, there are still a couple of open questions. Solving the conjecture on optimal couplings between two 
Poisson processes might help solve some of them. In [LT09 Last and Thorisson studied equivariant trans- 
ports between random measures in a rather general setting. In the recent article |LMT11| Last, Morters 
and Thorisson constructed an equivariant transport between two diffuse random measures to study un- 
biased shifts of Brownian motion. They also derive some moment estimates on the typical transport 
distance. Fair allocations have been studied and constructed, e.g. in [HP05 |HHP06, CPPRa l IMTllj 
and references therein. However, |HS10] is, to our knowledge, so far the only article studying couplings 
of two random measures under the additional requirement of being cost minimizing, that is optimal. 

1.1. Outline. In section [2] we introduce the setting and objects we work with. Section [3] is devoted to 
the proof of a key technical lemma, the existence and uniqueness of optimal semicouplings on bounded 
sets. In section [4] we prove Theorem |1.1| Theorem 1.2 and Theorem 1.3 are proved in section [5j In all 



these sections we always assume that the second marginal has intensity f3 < 1. In section [6] we treat the 
case of p > 1. In section [7] we state the estimates on the compound Poisson process. Finally in section [8] 
we show that W p defines a metric on V p and prove the stability result. 



2. Set-up 



In this section we will explain the general set-up, some basic concepts and derive the first result, a general 
existence result by a compactness argument. 
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2.1. The setting. From now on we will always assume to work in the following setting. (M,d,m) will 
denote a complete connected smooth non-compact Riemannian manifold with Riemannian distance d and 
Riemannian volume measure to. The Borel sets on M will be denoted by B(M). Given a map S and a 
measure p we denote the push forward of p under S by S*p, i.e. S*p(A) = p(S~ 1 (A)) for any Borel set A. 
Given any product X — H™ =1 Xi of measurable spaces, the projection onto the i-th space will be denoted 
by 7T.;. Given a set A C M its complement will be denoted by CA and the indicator function of A by I a- 

We will assume that there is a group G of isometries acting on M. For a set A C M we write r g A :— gA — 
{ga : a £ A}. For a point x E M its orbit under the group action of G is defined as Gx = {gx : g E G}. 
Its stabilizer is defined as G x = {g E G : gx — x} the elements of G that fix x. 

Definition 2.1 (Group action). Let G act on M. We say that the action is 

• properly discontinuous if for any x E M and any compact K C M gx E K for only finitely many 
gEG. 

• cocompact if M/G is compact in the quotient topology. 

• free if gx — x for one x E M implies g = id, that is the stabilizer for every point is trivial. 

We will assume that the group action is properly discontinuous, cocompact and free. By Theorem 3.5 in 
Bow06] this already implies that G is finitely generated and therefore countable. 

Definition 2.2 (Fundamental region). A measurable subset B C M is defined to be a fundamental 
region for G if 

i) U 9 eG9Bo = M 

ii) B n gB Q = for all id ± g E G. 

The family {gBo : g E G} is also called tessellation of M. 

There are many different choices of fundamental regions. We will choose a special one, namely a certain 
subset of the Dirichlet region with respect to some fixed point p. However, each fundamental region 
has the same volume and therefore defines a tiling of M in pieces of equal volume. Indeed, we have the 
following Lemma. 

Lemma 2.3. Let F\ and F 2 be two fundamental regions for G. Assume m(Fi) < 00. Then m(F±) = 
m(F 2 ). 

Proof. As F\ n gF 2 and F± D hF 2 are disjoint for g 7^ h by the defining property of fundamental regions 
we have 

m(Fi) = m (Fi n gF 2 ) = mig^Fi n F 2 ) = m{F 2 ). 

g&G 9 SG 

□ 

By scaling of the volume measure m we can assume that to (So) = 1- This assumption is just made to 
simplify some notations. 

As G is finitely generated, there are finitely many elements 01, . . . , a& E G such that every gEG can be 
written as a word in these letters and their inverses. The set S — {01, . . . , &&} is called a generating set. 
The generating set is not unique, e.g. Z is generated by {1} but also by {2,3}. We will fix one finite 
generating set for G. It does not matter which one as the results will be independent from the specific 
choice. 

Given the generating set S. We can construct a graph A = A(G, S) as follows. Put V(A) = G as the 
vertices. For each gEG and a E S we connect g and ag by a directed edge labeled with a. The same 
edge with opposite orientation is labeled by a -1 . This gives a regular graph of degree 215*1. We endow A 
with the word metric d& which coincides with the usual graph distance. 

Definition 2.4 (Cayley graph). If S is a generating set of G, then A(G,S) is called Cayley graph of G 
with respect to S. 

We denote the closed 2 r neighbourhood of the identity element in A by A r , that is A r = {gEG: 
^a(L 9) < 2 r }. The boundary of A r is defined as dA r = {h ^ A r : 3g E A r s.t. dA(h,g) = 1}. By B r 
we denote the range of the action of A r on the fundamental domain i?o, that is B r — UgeA 9^0- 

We will need to control the mass that is close to the boundary of B r , that is the growth of B r . In section 
[5j we will assume that 

urn r— = V g E G. 

r->oo |A r | 
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Several times we will use a rather simple but very powerful tool, the mass transport principle. It already 
appeared in the proof of Lemma |2.3| It is a kind of conservation of mass formula for invariant transports. 

Lemma 2.5 (mass transport principle). Let f : G x G — > R+ be a function which is invariant under the 
diagonal action of G, that is f[u,v) = f(gu,gv) for all g,u,v G G. Then we have 

veG veG 

Proof. 

^2f(u,v) = ^2f{u,gu) = ^2f{g- 1 u,u) = ^/(u,u). 

veG geG gGG veG 

□ 

For a more general version we refer to [BLPS99] and |LT09| . 

Recall the disintegration theorem for finite measures (e.g. see Theorem 5.1.3 in |AGS08j or 111-70 in 
[DM78p . 

Theorem 2.6 (Disintegration of measures). Let X, Y be Polish spaces, and let 7 be a finite Borel measure 
on X x Y . Denote by fi and v the marginals of 7 on the first and second factor respectively. Then, there 
exist two measurable families of probability measures ('y x )xex md ("fy) y eY such that 

l(dx,dy) = j x (dy)fi(dx) = j y (dx)v(dy). 

2.2. Couplings and Semicouplings. For each Polish space X (i.e. complete separable metric space) 
the set of Radon measures on X - equipped with its Borel er-ficld - will be denoted by M.(X). Given 
any ordered pair of Polish spaces X, Y and measures A G Ai(X),fi G A4(Y) we say that a measure 
q G A4(X x Y) is a semicoupling of A and fx, briefly q G Tl s (X,n), iff the (first and second, resp.) 
marginals satisfy 

(vri)*g < A, {^2)*q = H, 
that is, iff q(A x Y) < X(A) and q(X x B) = (j,(B) for all Borel sets Ac X,B cY. The semicoupling q 
is called coupling, briefly q £ II(A,/x), iff in addition 

(•7Tl)*<? = A. 

See also |FiglO| for the related concept of partial coupling. 

2.3. Random measures on M. We endow Ai(M) with the vague topology. The next Lemma summa- 
rizes some basic facts about vague topology (e.g. see [Kal97j or }Bau01| ) 

Lemma 2.7 (vague topology). Let X be a locally compact second countable Haussdorff space. Then, 

i) A4(X) is a Polish space in the vague topology. 

ii) A C M(X) is vaguely relatively compact iff sxvp fieA fJ-(f) < 00 for all f £ C C (X). 
hi) // /j, n A fi and B C X relatively compact with fj,(dB) = then ii n (B) ^{B). 

The action of G on M induces an action of G on Ai(M x . . . x M) by push forward with the map r g : 

(T g ) m X(A u ...,A k ) = A(( ff - 1 (A 1 ),..., ff - 1 (A fe )) yA u ...A k eB{M),ke'N. 

A random measure on M is a random variable A* (the notation with the "•" is intended to make it easier 
to distinguish random and non-random measures) modeled on some probability space (fl, 21, P) taking 
values in M.(M). It can also be regarded as a kernel from f2 to M. Therefore, we write either X bJ (A) 
or X(uj,A) depending on which property we want to stress. For convenience, we will assume that Q is 
a compact metric space and 21 its completed Borel field. These technical assumptions are only made to 
simplify the presentation. 

A point process is a random measure \x* taking values in the (vaguely closed) subset of all locally finite 
counting measures on M. It is called simple iff ^({x}) G {0, 1} for every x G M and a.e. lo G O. We 
call a random measure A* absolutely continuous iff it is absolutely continuous to the volume measure m 
on M for a.e. u> G O. It is called diffusive iff it has no atoms almost surely. The intensity measure of a 
random measure A* is a measure on M defined by A 1— > E[A*(^4)]. 

The class of all relatively compact sets in B(M) will be denoted by B. For a random measure A* its class 
of stochastic continuity sets is defined by B\> — {A G B : X'(dA) = a.s.}. Convergence in distribution 
and tightness in Ai(M) can be characterized by 
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Lemma 2.8 (tightness of random measures). Let A*, A*, . . . be random measures on M. Then the sequence 
(A*) n£ N is relatively compact in distribution iff (A* (A)) n ^ is tight in R + for every A G B. 

Theorem 2.9 (convergence of random measures). Let A*, A*, A*,... be random measures on M. Then, 
these conditions are equivalent: 

i) A* 4 A' 

ii) A^(/)4a'(/) for all f eC c (M) 

hi) (A* (Ax), . . . , X' n (A k )) 4 (A'(Ai), . . . , X'(A k )) for all A u . . . ,A k G B x .,k G N. 

If A* is a simple point process or a diffusive random measure, it is also equivalent that 
iv) A* (A) 4 A* (A) for all A G B\>. 

For the proof of these statements we refer to Lemma 14.15 and Theorem 14.16 of |Kal97j . 

Just as in Lemma 11.1. II of |DVJ07j we can derive the following result on continuity sets of a random 

measure A*: 

Lemma 2.10. Let A* be a random measure on M, A G B(M) be bounded and (A) r be the r-neighbourhood 
of A in M. Then for all but a countable set of r G R+ we have (A) r G B\> . 

A random measure A* : O —> A4(M) is called G-invariant or just invariant if the distribution of A* is 
invariant under the action of G, that is, iff 

(r 9 )*A* 4' A* 

for all g G G. A random measure q* : O — > A4(M x M) is called invariant if its distribution is invariant 
under the diagonal action of G. 

If (0,21) admits a measurable flow 9 g : f2 0, g G G, that is a measurable mapping (ui,g) >-> O g uj with 
6q the identity on ft and 

6 g o 9 h = 6 gh , g, he G, 
then a random measure A* : O — > A / t(M) is called G-equivariant or just equivariant iff 

A(0 9 w,c/A) = X(u,A), 

for all .g G G,u G G B(M). We can think of \(6 g oj,-) as \(w, ■) shifted by g. Indeed, let SOt be 

the cylindrical <r— algebra generated by the evaluation functionals A i— > ^(A), A G B(M),[i e M. As in 
example 2.1 of |LT09| . consider the measurable space (A^,9Jt) and define for /j e M,g G G the measure 
9 g fi(A) = ji(g~ 1 A). Then, {9 g ,g G G} is a measurable flow and the identity is an equivariant measure. 
A random measure q' : fl — > M.[M x M) is called equivariant iff 

q 6 ° u (gA,gB) = q u (A,B), 

for all g G G, w G ft, A, B G B(M). 

Example 2.11. Let q* be an equivariant random measure on M x M given by = (id, T")*A" for 
some measurable map T* and some equivariant random measure A*. The equivariance condition 



l B (y)5 T e g ^ gx) (d(gy))\°^(dx) = q°>»(gA,gB) = (f(A,B) - / l B (y)6 T « {x) (dy)\"(dx), 

A J A 

translates into an equivariance condition for the transport maps: 

T e °"(gx)=gT»(x). 
A probability measure P is called stationary iff 

¥o9 g =¥ 

for all g G G. Given a measure space (0,21) with a measurable flow {O g ) g eG an d a stationary probability 
measure P any equivariant measure is automatically invariant. The advantage of this definition is that the 
sum of equivariant measures is again equivariant, and therefore also invariant. The sum of two invariant 
random measures does not have to be invariant (see Remark |2.19 1. 



We say that a random measure A* has subunit intensity iff E[A* (A)] < m(A) for all A G B(M). If equality 
holds in the last statement we say that the random measure has unit intensity. An invariant random 
measure has subunit (or unit) intensity iff its intensity 

(3 = E[X'(B )} 
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is < 1 (or = 1 resp.). Given a random measure, the measure (X*P)(dy, doj) := X UJ (dy)¥(duj) on M X fi is 
called Campbell measure of the random measure A*. 



Example 2.12. i) The Poisson point process with intensity measure m. It is characterized by 

— for each Borel set A C M of finite volume the random variable uj h-> ^(A) is Poisson 
distributed with parameter m(A) and 

— for disjoint sets Aj,...^ C M the random variables (i u (Ax), . . . , ^(A/,) are independent. 
It can be written as 

with some countable set C M without accumulation points, 
ii) The compound Poisson process is a Poisson process with random weights instead of unit weights. 
It is compounded with another distribution giving the weights of the different atoms. It can be 
written as 

E x & 

££H(u;) 

for some iid sequence (-2Q)geE(uj) independent of the Poisson point process. For example one 
could take X^ to be a Poisson random variable or an exponentially distributed random variable. 
If has distribution 7 we say fi* is a 7— compound Poisson process. 

From now on we will always assume that we are given two equivariant random measures A* and /i* 
modeled on some probability space (fi, 21, P) admitting a measurable flow {0 g ) ge c such that P is stationary. 
We will assume that fi is a compact metric space. Moreover, we will assume that A* is absolutely 
continuous and A* and \i* are almost surely not the zero measure. Note that the invariance implies that 
H U (M) = \"(M) = 00 for almost every w (e.g. see Proposition 12.1.VI in |DVJ07j ). 

2.4. Semicouplings of A* and fi*. A semicoupling of the random measures A* and fi' is a measurable 
map q* : fi -> M(M x M) s.t. for P-a.e. u e fi 

q M is a semicoupling of A" and /i". 

Its Campbell measure is given by Q — q'F € A4(M X M X fi). Q is a semicoupling between the Campbell 
measures A'P and /i'P in the sense that 

Q(M x • x ■) = /i'P and Q(< x M x •) < A'P. 

Q could also be regarded as semicoupling between A'P and /i'P on M x fi x M x fi which is concentrated 
on the diagonal of fi x fi. It could be interesting to relax this last condition on Q and allow different 
couplings of the randomness. However, we will not do so and only consider semicouplings of A'P and 
/i'P that are concentrated on the diagonal of fi x fi. We will always identify these semicouplings with 
measures on M x M X fi. 



Given such a semicoupling Q G M.{M x M x fi) we can disintegrate (see Theorem 2.6) Q to get a 
measurable map q' : fi — > Ai(M x M) which is a semicoupling of A* and /i*. 

According to this one-to-one correspondence between q* — semicoupling of A* and /i* — and Q = q'F 
— semicoupling of A'P and /i'P — we will freely switch between them. And quite often, we will simply 
speak of semicouplings of A* and fi' . 

We denote the set of all semicouplings between A* and /z* by II s (A*,/x*). The set of all equivariant 
semicouplings between A* and fi' will be denoted by II es (A*,/i*). 

A factor of some random variable X is a random variable Y which is measurable with respect to &(X). 
This is equivalent to the existence of a deterministic function / with Y = f(X). In other words, a factor 
is a rule such that given X we can construct Y. A factor semicoupling is a semicoupling of A* and /i* 
which is a factor of A* and \i* . 

2.5. The Monge-Kantorovich problem. Let A, /i be two probability measures on M. Moreover, let 
a cost function c : M x M R be given. The Monge-Kantorovich problem is to find a minimizer of 

c(x,y) q(dx,dy) 

MxM 

among all couplings q of A and /i. A minimizing coupling is called optimal coupling. If the optimal 
coupling q is induced by a transportation map, i.e. q = (id,T)*\, we say that q is a solution to the 
Monge problem. There are rather general existence and uniqueness results for optimal couplings. We 
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always assume that the cost function c(x,y) = d{d(x,y)) is such that there is a unique solution to the 
Monge problem between A and \i whenever A <C m. For conditions on d such that this assumption is 
satisfied we refer to section [3] 

It can be shown that any optimal coupling is concentrated on a c— cyclical monotone set. A set A C XxX 
is called c— cyclical monotone if for all n £ N and (xi,yi)f =1 € A n it holds that 

n n 

^2,c{x l ,y l ) < ^2c(xi,y l+1 ), 

i=l i—1 

where y\ = y n +i- If the cost function is reasonably well behaved (continuous is more than sufficient, 
see BGMS09 ), also the reverse direction holds. Any coupling which is concentrated on a c— cyclical 
monotone plan is optimal. For further details and applications of mass transport theory we refer to 
[RR98llViItklVilQ9] 

2.6. Cost functionals. Throughout this article, ■d will be a strictly increasing, continuous function from 
R + to M + with $(Q) = and lim l d(r) — oo. Given a scale function $ as above we define the cost function 

r— >oo 

c(x,y) =d(d(x,y)) 

on M x M, the cost functional 

Cost(g) = / c(x, y) q(dx, dy) 

JMxM 

on M (M x M) and the mean cost functional 

£ost(Q) — / c(x,y) Q(dx,dy,duj) 

JMxMxU 

on M(M x M x fi). 

We have the following basic result on existence and uniqueness of optimal semicouplings the proof of 
which is deferred to section [3] The first part of the theorem, the existence and uniqueness of an optimal 
semicoupling, is very much in the spirit of an analogous result by Figalli FiglO| on existence and (if 
enough mass is transported) uniqueness of an optimal partial coupling. However, in our case the second 
marginal is arbitrary whereas in |FiglO| it is absolutely continuous. 

Theorem 2.13. (i) For each bounded Borel set A C M there exists a unique semicoupling Qa of A*P 
and (l,4/i*)P which minimizes the mean cost functional £ost(.). 

(ii) The measure Qa can be disintegrated as Q A(dx , dy , duS) :— q^dx, dy) P(dw) where for ¥ '-a. e. uj the 
measure q 1 ^ is the unique minimizer of the cost functional Cost(.) among the semicouplings of A" and 
1 >/' • 

(Hi) £ost(Q A ) = J Cost(g^)P(dw). 

For a bounded Borel set A C M , the transportation cost on A is given by the random variable : f2 — > 
[0, oo] as 

C a (oj) := Cost(q^) = inf{Cost(g w ) : q u semicoupling of A" and Ia^}- 
Lemma 2.14. (i) If A\, . . . ,A n are disjoint then Vcj € fi 

4=1 1 i = l 

(ii) If A\ = gA2 for some g € G, then Ca x and Ca 2 o-re identically distributed. 

Proof. Property (ii) follows directly from the joint invariance of A* and /x*. The intuitive argument for 
(i) is, that minimizing the cost on IJ^ Ai is more restrictive than doing it separately on each of the Ai. 
The more detailed argument is the following. Given any semicoupling q w of A" and ly . a^ then for 
each % the measure qf := 1 MxA,- q u is a semicoupling of A w and lytiM w - Choosing q" as the minimizer of 
C » (u>) yields 

i = l 

C UiA M = Cost(^) =^Cost(gf) > ^C^Cw). 

i i 

□ 
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2.7. Standard tessellations. In this section, we construct the fundamental region B and thereby a 
tessellation or a tiling of M. We will call this tessellation a standard tessellation. The specific choice of 
fundamental domain is not really important for us. However, we will choose one to fix ideas. 

We now define the Dirichlet region. To this end let p £ M be arbitrary. Due to the assumption of 
freeness, the stabilizer of p is trivial. Construct the Voronoi tessellation with respect to Gp, the orbit of 
p. The cell containing p is the Dirichlet region. 

Definition 2.15 (Dirichlet region). Let p G M be arbitrary. The Dirichlet region of G centered at p is 
defined by 

D P (G) = {xeM : d{x,p) <d(x,gp)Vg G G}. 
From now on we will fix p and write for simplicity of notation D = D p {G). We want to construct a 

o 

fundamental domain from D. For every x G D we have d(x,p) < d(gx,p) for every id ^= g G G, that is 
\Gx fl D\ = 1, where \H\ denotes the cardinality of H. However, ii x £ 3D we have x G D n gD ^ for 
some g G G. This implies that \Gx H D\ > 2. Yet, for the fundamental region, Bq, we need exactly one 
representative from every orbit. Hence, we need to chose from any orbit Gx intersecting the boundary 
of D exactly one representative z G Gx n 3D. Let V be a measurable selection of these and finally define 

o 

-Bo = D U V. By definition, Bq is a fundamental region. Such a measurable selection exists by Theorem 
17 and the following Corollary in |Del75j . 

Example 2.16. Considering M. d with group action translations by Z d a choice for the fundamental region 
would be Bq — [0, l) d . If we consider M = H 2 the two dimensional hyperbolic space we can take for G 
a Fuchsian group acting cocompactly and freely, that is, with no elliptic elements. Then, the closure of 
the Dirichlet region becomes a hyperbolic polygon (see |Kat92j ). 

2.8. Optimality. The standard notion of optimality - minimizers of Cost or Cost - is not well adapted 
to our setting. For example for any semicoupling q* between the Lebesgue measure and a Poisson point 
process of intensity /3 < 1 we have £ost(q*) = oo. Hence, we need to introduce a different notion which 
we explain in this section. 

The collection of admissible sets is defined as 

Adm(M) = {B G B(M) : 31 C G, 1 < |/| < oo, F fundamental region : B = \J gF}. 

For a semicoupling q' between A* and fj,* the mean transportation cost of q* is defined by 
£(0 := sup — 



c{x,y) q*(dx,dy) 



MxB 

Definition 2.17. A semicoupling q* between A* and fi* is called 

i) asymptotically optimal iff 

£(<?•)= inf £(?•) =: c^. 

ii) optimal iff q* is equivariant and asymptotically optimal. 
We will also use several times the quantity 

inf £(<?•) =: Coo . 

Obviously < c £y00 . 

Note that the set of optimal semicouplings is convex. This will be useful for the proof of uniqueness. 
Remark 2.18. Equivariant semicouplings q* are invariant. Hence, they are asymptotically optimal iff 

£(q')=E / c{x,y)q'{dx,dy) = c ejOC .. 

J My. B 

Because of the invariance, the supremum does not play any role. Moreover, for two different fundamental 
regions Bq and Bq define 

f(g,h) = E[Cost{l Mx(gBonhSo) q')}. 

Then, for k G G and equivariant q* we have f(g, h) — f(kg, kh). Hence, we can apply the mass transport 
principle to get 



E 



c(x,y)q'(dx,dy) 



M x S 



^2 f(id,h) = f(g,id) / c{x,y)q'(dx,dy) 

heG g£G UMxBo 
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Thus, the specific choice of fundamental region is not important for the cost functional (£(•) if we restrict 
to equi variant semicouplings. 

Remark 2.19. The notion of optimality explains why we restrict to stationary probability measures and 
equivariant random measures. If A* and p* are just invariant, there does not have to be any invariant 
semicoupling between them. Indeed, take A* a Poisson point process of unit intensity in M. d . It can be 
written as p, u — Yl^estu) % Define A" := X^eH(u) to ^ e * ne Poisson process that we get if we reflect 
the first one at the origin. Then A* and p* are invariant but not jointly invariant, e.g. consider the set 
[0, l) d x [— 1,0) , and not both of them can be equivariant. 

2.9. An abstract existence result. Given that the mean transportation cost is finite the existence of 
an optimal semicoupling can be shown by an abstract compactness result. A similar reasoning is used to 
prove Corollary 11 in [Hol09 . 

Proposition 2.20. Let A* and p,* be two equivariant random measures on M with intensities 1 and 
j3 < 1 respectively. Assume that int^gri^A*,^*) £(<Z*) — c e,oo < oo, then there exists some equivariant 
semicoupling q* between A* and p* with £(q°) = c e;00 . 

Proof. As c e .oo < oo there is a sequence g* £ ir s (A* , p') such that (£(<?*) — c n \ c ej00 . Moreover, we can 
assume that the transportation cost is uniformly bounded by c„ < 2c e ,oo =: c for all n. We claim that 
there is a subsequence (q' lk )k<£N of (<?*)neN converging to some q* 6 IIj s (A*,/x*) with <t(q') — c e ,oo- We 
prove this in four steps: 



i) The functional £(•) is lower semicontinuous: 

It is sufficient to prove that the functional Cost(-) is lower semicontinuous. Let (p n )neN De an Y sequence of 
couplings between finite measures converging to some measure p in the vague topology. If Cost(p„) = oo 
for almost all n we are done. Hence, we can assume, that the transportation cost are bounded. Let [Bq) t 
denote the r-neighbourhood of B . For k 6 K let <^ : M x M — > [0,1] be nice cut off functions with 
<f) k (x, y) = l on (B ) k x (B ) k and (f> k {x, y) = if x G C((B )fc+i) or y e C((B )fe+i). Then, we have using 
continuity of the cost function c(x,y) and by the definition of vague convergence 



liminf Cost(p„) = liminf / c(x,y)p n (dx,dy) 

n^co ' n^oo J MxM 



= lim inf sup / 4>k(x,y) c(x,y)p n (dx, dy) 

n^roo ken J My. M 

> sup lim inf / (f> k (x,y) c(x,y)p n (dx,dy) 
k n ^°° JMxM 



= sup/ <f>k(x,y) c(x,y)p(dx,dy) = Costco). 

fc JMxM 

Applying this to lAfxs 9* shows the lower semicontinuity of £(•). 

ii) The sequence (q') ne n is tight in M(M x M x f2): 

Put / € C C (M x M x CI). According to Lemma 2.7 we have to show sup neN q*P(/) < Mt < oo for some 
constant Mf. To this end let A C M compact be such that supp(/) C A x M x f2. We estimate 

f f(x,y,uj)q%(dx,dy)¥(duj) < \\f\\oo A*P(^4 x Q) 

J MxMxfl 

< H/lloo m(A) =: Mf. 

Hence, there is some measure q* and a subsequence q* nh with q* fc — > q* in vague topology on A4(M x 
M x fi). By lower semicontinuity, we have £(q') < liminf ) = c eiOC . Now we have a candidate. We 
still need to show that it is admissible. 

iii) q' is equivariant: 

Take any continuous compactly supported / g C C (M x M x fl). By definition of vague convergence 



f(x,y,Lj)q^ h (dx,dy)F(du)^ J f(x,y,u)q"(dx,dy)P(du). 
As all the g* fc are equivariant, we have for any g € G 

f(x,y,u)£ h (dx,dy)F(u) = f f(g- 1 x,g- 1 y,e g u;)q e n f(dx,dy)P(du) 

f(g- 1 x,g~ 1 y,e g oj)q e ^(dx,dy)¥(duj) 
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f(g- L x, g- L y, 6 g Lu)q ti ^ (dx, dy)P(duj). 



Putting this together, we have for any g £ G 

f{x,y,uj)q UJ (dx,dy)V{dLo) = 
Hence, q* is equi variant. 

iv) q* is a semicoupling of A* and fj,' : 

Fix h £ C C (M x 0). Put A C M compact such that supp(fr) C A x SI and A £ Adm(M). Denote the 
R— neighbourhood of A by A R . By the uniform bound on transportation cost we have 

(2-1) q^>(C(A R ),A,Sl)<m(A)^ 

uniformly in n. Let f R : M — ¥ [0, 1] be a continuous compactly supported function such that f R {x) = 1 
for x £ An and f R {x) = for x £ CA R+ i. As q*P is a semicoupling of A* and [i° we have due to monotone 
convergence 



%,wK(dy)P(dw) 



A/xSl 



%,u;)#(da;,dj/)P(<M 



MxMxfi 



lim 



f R (x)h(y,w)q%(dx,dy)W(duj). 



M xJlfxfl 



Because of the uniform bound (2.1 1 we have 



/i(x,w)/x"(dx)P(dw) 



m x n 



< m(A) 



c-\\h\\ 



f R (x)h(y,u;)q% k (dx,dy)F(du>) 

I MxMxQ 

Taking first the limit of — > oo and then the limit of R —¥ oo we conclude using vague convergence and 
monotone convergence that 



= lim lim 

R — ^oo k — s-oo 



h(y,uj)^(dy)F(du)- 



lim 

R— >oo 



MxQ 

h(y,oj)^(dy)P(du;) 

M x Q 

h(y,u;) f x u (dy)¥(du)- 



f R (x)h(y,oj)q^ k (dx,dy)P(du J ) 

MxMxQ 

f R (x)h(y,u;)q u (dx,dy)P((kj) 



MxQ 



M x M x n 

h(y,uj)q u} (dx,dy)¥(duj) 



M x M x n 



This shows that the second marginal equals /z*. For the first marginal we have for any k £ C C {M x SI) 



MxQ. 



k(x,u)q% k (dx,dy)V(<k;) < 



k(x,oj)X u (dx)¥(du;). 



MxQ. 



MxQ 



k(x,uj)X ul (dx)V{duj) 



MxQ 



In particular, using the function f R from above we have, 

f R (y) k(x,Lj)q% h (dx,dy)F(dw) < 
Taking the limit — > oo yields by vague convergence 

f R (y) k(x,oj)q^(dx,dy)¥(du J ) < 

MxQ 

Finally taking the supremum over R shows that q* is indeed a semicoupling of A* and /i*. 



k(x,u})\ u (dx)W>(du). 



m x n 



□ 



Remark 2.21. i) This coupling need not be a factor coupling. We do not know if it is in general 

true or not that = c e ,oo, that is, if minimizing the functional <£.(■) over all semicouplings is 
the same as minimizing over all equivariant semicouplings. However, in the case that the balls 
A r cG are F0lner sets, we can show equality (see Corollary [5^5] and Remark 5.6 1, 
ii) The same proof shows the existence of optimal semicouplings between A* and /x* with intensities 
1 and (3 > 1 respectively. In this case the "semi" is on the side of fj,' (see also section [6]). 

Lemma 2.22. Let q* be an invariant semicoupling of two random measures A* and \x* with intensities 
1 and (3 < 1 respectively. Then, q* is a coupling iff (3 = 1 . 

Proof. This is another application of the mass transport principle. Let Bq be a fundamental region and 
define f(g, h) — E[q*(gBo, hBo)]. By invariance of q', we have f(g, h) — f(kg, kh) for any k £ G. Hence, 
we get 

1 = E[A'(S )] > E[q'(B ,M)] = £ f(id,g) = £ f(h, id) = E[q'(M,B )] = [3. 

g£G h£G 
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We have equality iff j3 = 1. By definition of semicoupling, we also have q^(A,M) < X u (A) for any 
A C M. Hence, in the case of equality we must have q u (A,M) = A^(^4) for P— almost all u). □ 

Remark 2.23. The remark above applies again. Considering the case of intensity f3 > 1 gives that q' is 
a coupling iff j3 = 1. 

2.10. Assumptions. Let us summarize the setting and assumptions we work with in the rest of the 
article. 

• M will be a smooth connected non-compact Riemannian manifold with Riemannian volume mea- 
sure m, such that there is a group G of isometries acting properly discontinuously, cocompactly 
and freely on M. 

• Bq will denote the chosen fundamental region. 

• c(x,y) = $(d(x,y)) with : R + — > R + such that {)(0) = and \im r ^ OD , d(r) — oo. Given two 
compactly supported probability measures on M, A <C m and p arbitrary, we will assume that 
the optimal transportation problem admits a unique solution which is induced by a measurable 
map T, i.e. q = (id, T)*A. 

• (0,21, P) will be a probability space admitting a measurable flow (d g ) g€ c- P is assumed to be 
stationary and f2 is assumed to be a compact metric space. 

• A* and p* will be equivariant measure of intensities one respectively f3 E (0,oo). Moreover, we 
assume that A* is absolutely continuous. 



3. Optimal Semicouplings on bounded sets 

The goal of this section is to prove Theorem |2.13| the crucial existence and uniqueness result for optimal 
semicouplings between A* and p* restricted to a bounded set. The strategy will be to first prove existence 
and uniqueness of optimal semicouplings q — q u for deterministic measures A = A" and p — p? . Secondly, 



we will show that the map uj i— > q" is measurable, which will allow us to deduce Theorem 2.13 



Optimal semicouplings are solutions of a twofold optimization problem: the optimal choice of a density 
p < 1 of the first marginal A and subsequently the optimal choice of a coupling between pX and p. This 
twofold optimization problem can also be interpreted as a transport problem with free boundary values. 
Throughout this section, we fix the cost function c(x,y) — {>(d(x,y)) with ■& - as before - being a 
strictly increasing, continuous function from R + to R + with $(0) = and lim §(r) = oo. As already 

r— >oo 

mentioned, we additionally assume that the optimal transportation problem between two compactly 
supported probability measures A and p such that A <§; to has a unique solution given by a transportation 
map, e.g. the optimal coupling is given by q = (id, T)*A. There are very general results on the uniqueness 
of the solution to the Monge problem for which we refer to chapters 9 and 10 of |Vil09j . To be more 
concrete we state a uniqueness result for compact manifolds due to McCann jMcCOlj and an uniqueness 
result by Huesmann and Sturm in the simple but for us very interesting case that the measure p is 
discrete. 

Theorem 3.1 (McCann). Let N be a compact manifold, A <C m and p be probability measures and 
c(x,y) = j^ x,v ^ r(s)ds with r : M.+ — > R continuously increasing and r(0) = 0. Then, there is a 
measurable map T : M M U {9} such that the unique optimal coupling between A and p is given by 
q = (id, T)*A. 

The "cemetery" 3 in the statement is not really important. This is the place where all points outside of 
the support of A are sent. We just include it to make some notations easier. 

If we assume p to be discrete, Lemma 6.1 in [HS10 shows that we can actually take $ to be any continuous 
strictly increasing function. 

Lemma 3.2. Given a finite set 5 = {£i, . . . C M, positive numbers ((Xi)i<j<fc summing to one and 
a probability density p 6 L 1 (M, m). Consider the cost function c(x,y) = i)(d(x,y)) for some continuous 
strictly increasing function $ : M + — > K + such that $(0) = and lim^oo -d(r) — oo. If dim(M) = 1 we 
exclude the case $(r) = r. 

i) There exists a unique coupling q of p ■ m and a = X^=i fli<% which minimizes the cost functional 
Cost(-). 

ii) There exists a (m-a.e. unique) map T : {p > 0} — > 3 with T*(p ■ m) — a which minimizes 
J c(x, T(x))p(x) m(dx). 

Hi) There exists a (m-a.e. unique) map T : {p > 0} — > E with T„(p- to) = a which is c-monotone (in the 
sense that the closure of {(x,T(x)) : p(x) > 0} is a c-cyclically monotone set). 
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iv) The minimizers in (i), (ii) and (Hi) are related by q — (Id,T)*(p ■ m) or, in other words, 

q(dx,dy) = 5 T ( x ){dy) p(x)m(dx). 

Remark 3.3. In the case that dim(M) — 1 and cost function c(x,y) = d(x,y) the optimal coupling 
between an absolutely continuous measure and a discrete measure need not be unique. In higher dimen- 
sions this is the case, as we get strict inequalities in the triangle inequalities. A counterexample for one 
dimension is the following. Take A to be the Lebesgue measure on [0, 1] and put p = |i5o + |^i/i6- Then, 
for any a e [1/16, 1/3] 

q a (dx,dy) = l [0 ^(x)5 (dy)X(dx) + l[ a ,2/3+a)(x)5 1/1G (dy)\(dx) + l[ a +2/3,i]{x)S (dy)X(dx) 
is an optimal coupling of A and p with Cost((7 a ) = 11/24. 

Remark 3.4. In the case of cost function c(x,y) — -d p (x,y) the optimal transportation map is given 

by 

TM - exp(<*M)|||! 

for functions = ~^ i d p (z,^ i ) + with constants bi and j such that &j(x) = maxKKi: ^i(x) (e.g. 

see jMcCOlp . 

Given two deterministic measures A = f-m for some compactly supported density / (in particular A <C m) 
and an arbitrary finite measure p with supp(/i) C A for some compact set A such that p(M) < X(M) < oo. 
We are looking for minimizers of 

Cost(g) = J c(x, y)q(dx, dy) 

under all semicouplings q of A and p. The key step is a nice observation by Figalli, namely Proposition 
2.4 in |FiglO| . The version we state here is adapted to our setting. 

Proposition 3.5 (Figalli). Let q be a Cost minimizing semicoupling between A and p. Write f q -m = 
(tti)*9- Consider the Monge-Kantorovich problem: 

minimize C '(7) = / c(x,y) n f(dx,dy) 

J MxM 

among all 7 which have A and fj, + (/ — f q ) ■ m as first and second marginals, respectively. Then, the 
unique minimizer is given by 

q + (id x id)*(f - f q ) ■ m. 

This allows us to show that all minimizers of Cost are concentrated on the same graph which also gives 
us uniqueness: 

Proposition 3.6. There is a unique Cost minimizing semicoupling between A and fx. It is given as 
q = (id, T)*(p ■ A) for some measurable map T:M-}MU {3} and density p. 

Proof, (i) The functional Cost(-) is lower semicontinuous on A4(AI x M) wrt weak convergence of mea- 
sures. Indeed, take a sequence of measures (q n )n&i converging weakly to some q. Then we have by 
continuity of the cost function c(-, •): 

c(x,y) q(dx,dy) — sup / c(x , y) A k q(dx , dy) 

keNJ 



= sup lim / c(x , y) A k q n (dx , dy) 

keN n^ocJ 

< liminf / c(x,y) q n (dx,dy). 

n— > 00 / 

(ii) Let O denote the set of all semicouplings of A and p and 0\ denote the set of all semicouplings q 
satisfying Cost(g) < 2 inf^go Cost(q) =: 2c. Then 0\ is relatively compact wrt weak topology. Indeed, 
q(M x CA) = for all q G Q x and 

q(C(A r ) x A) < J-j ■ Cost(g) < 

for each r > where A r denotes the closed r-neighborhood of A in M. Thus, for any e > there exists 
a compact set K = A r x A in M x M such that g(C/v) < e uniformly in q € 0\. 

(iii) The set O is closed wrt weak topology. Indeed, if q n — > q then (iTi) t q n — > (TTi)*q and (7r 2 )*(j , n — > 
(iT2)*Q- Hence 0\ is compact and Cost attains its minimum on O. Let q denote one such minimizer. Its 
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first marginal is absolutely continuous to m. By Theorem |3 . 1 1 and Lcmma |3.2[ there is a measurable map 
T : M ->■ M U {9} and densities /„ f q such that q = {id, !%(/, ■ A) = {id, T)„{f q ■ to). 



(iv) Given a minimizer of Cost, say q. By Proposition 3.5 q := q + {id,id) Sf {f — f q ) ■ m solves 

minC(7) = J c{x,y)j{dx,dy) 
under all 7 which have A and p + (/ — f q )m as first respectively second marginals, where f q -m = {"Ki)*q 



as above. By Theorem 3.1 and Lemma |3.2| there is a measurable map S such that q = {id, S%A. That 
is, q and in particular q are concentrated on the graph of S. By definition q = q + {id,id)*(f — f q ) ■ m 
and, therefore, we must have S{x) = x on {/ > f q }. 

(v) This finally allows us to deduce uniqueness. By the previous step, we know that any convex combi- 
nation of optimal semicouplings is concentrated on a graph. This implies that all optimal semicouplings 
are concentrated on the same graph. Moreover, Proposition |3.5| implies that if we do not transport all 
the A mass in one point we leave it where it is. Hence, all optimal semicouplings choose the same density 
p of A and therefore coincide. 

Assume there are two optimal semicouplings q± and q%. Then q 3 :— \{q\ + q 2 ) is optimal as well. By 
the previous step for any i £ {1,2,3}, we get maps Si such that is concentrated on the graph of Si- 
Moreover, we have 63(2) = x on the set {/ > f q3 } = {/ > / 91 }U{/ > f q2 }, where again f qi -to = (ni)*qi. 
As q 3 is concentrated on the graph of S3, q± and q 2 must be concentrated on the same graph. Hence, 
we have £3 = Si on {f q . > 0} for i = 1,2. We also know from the previous step that Si{x) = x on 
{/ > /»} C {/ > f q3 }. This gives, that S 3 = Si = S 2 on {/ > 0}. 

We still need to show that {f qi > 0} = {f q2 > 0}. Put A\ := {f qi > f q2 } and A 2 := {f q , 2 > f qi } and 
assume m(Ai) > 0. As A± C {/ > f q2 } we know that S 3 {x) = x on A± and similarly Ss(x) = x on A 2 . 
Now consider 

A := S^iAx) = (An{f qi =f q2 })U(AnAi)U(AnA 2 ). 
As S 3 (A 2 ) C A 2 and ^10/12 = we have A n A2 = 0. Therefore, we can conclude 
/*(i4i) = {S 3 )*f qi m{Ai) = f gi m(A 1 ) + f gi m{An{f qi = f qs }) 

> f q MM) + f q2 m{A n {f qi = f q2 }) 

= {S 3 )*f q2 m{A 1 ) = 

which is a contradiction, proving q± = q 2 . □ 

Remark 3.7. Let q = {id,T)*{p\) be the optimal semicoupling of A and p. If p happens to be discrete, 
we have p{x) £ {0, 1} m almost everywhere. Indeed, assume the contrary. Then, there is £ € supp(/x) 
such that on U := T _1 (£) we have p £ (0,1) on some set of positive A measure. Let R be such that 
\(UDB(£,R)) = p{{£}), where B{£,R) denotes the ball of radius R around £. Put V = f/n and 

q{dx,dy) = q{dx,dy) — \u{x)p{x)5^{dy)\{dx) + \v{x)5^{dy)\{dx). 

This means, we take the same transportation map, but use the A mass more efficiently q leaves some A 
mass far out and instead uses the same amount of A mass which is closer to the target £. By construction, 
we have Cost(g) > Cost(q) contradicting optimality of q. 

We showed the existence and uniqueness of optimal semicouplings between deterministic measures. The 



next step in the proof of Theorem 2.13 is to show the measurability of the mapping a; H> 3?(A W , 1a^ w ) = Qa 
the unique optimal semicoupling between A w and IaP^ ■ The mapping uj M- (A", IaP^) is measurable by 
definition. Hence, we have to show that (A w , IaP^) $(A", ^AP U ) is measurable. We will show a bit 
more, namely that this mapping is actually continuous. We start with a simple but important observation 
about optimal semicouplings. 

Denote the one-point compactification of M by M U {9} and let #(r) be such that it is equal to i?(r) 
on a very large box, say [0,K] and then tends continuously to zero such that c(x, 3) = ${d{x, 9)) = 
lim,.^oo #(f") = for any x £ M. By a slight abuse of notation, we also write 9 : M — > {9} for the map 
x 1 — ^ 9. 

Lemma 3.8. Let two measures A and p on M be given such that 00 > \{M) = N > p{M) = a 
and assume there is a ball B{x,K/2) such that supp{\) , supp{p) C B{x, K/2). Then, q is an optimal 
semicoupling between A and p wrt to the cost function c{-, •) iff q = q + {id, 9)*(1 — f q ) ■ A is an optimal 
coupling between A and p = p + {N — a)5g wrt the cost function c{-, •), where (fti)*q = f q X. 

Proof. Let q be any semicoupling between A and p. Then q = q + {id, 9)*(1 — f q ) ■ A defines a coupling 
between A and p. Moreover, the transportation cost of the semicoupling and the one of the coupling are 
exactly the same, that is Cost(g) = Cost(q). Hence, q is optimal iff q is optimal. □ 
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This allows to deduce the continuity of $ from the classical theory of optimal transportation. 

Lemma 3.9. Given a sequence of measures (\ n ) n &i converging vaguely to some X, all absolutely contin- 
uous to m with \ n {M) = X(M) = oo. Moreover, let (fx n )n£N be a sequence of finite measures converging 
weakly to some finite measure fx, all concentrated on some bounded set A C M. Let q n be the optimal 
semicoupling between A„ and fx n and q be the optimal semicoupling between X and fx. Then, q n converges 
weakly to q. In particular, the map (A,//) t— > $(A,/i) = q is continuous. 

Proof, i) As (fi n )neN converge to fx and fx is finite, we can assume that sup ra fx n (M), fi(M) < a < oo. 
As A n and A have infinite mass for any x G M and k G K there is a radius R{x, k) < oo such that 
X(B(x, R(x, k))) > k, where B(x, R) denotes the closed ball around x of radius R. Fix an arbitrary x € A 
and set R\ = R(x, a) + diam(A) and R2 = R(x, 2a) + diam(A). Because (A n )„ e N converge to A we can 
assume that for any n X n (B(x, R2)) > X(B(x, J?i) = N > a. 

ii) Optimality of q n and q implies that supp(g rl ) C B(x,R 2 ) x A and supp(g) C B(x,Ri) x A. Because 
otherwise there is still some mass lying closer to the target than the mass which is transported into the 
target. For any n let r n < R2 be such that X n (B(x,r n )) = N. Such choices exist as X n <C m for all n. 
Then, we even know that supp((7„) C B(x,r n ) x A. Set A„ = ls( X)r „)A„ and A = 1b(x,Ri)^- Then the 
optimal semicoupling between A ra and fi n is the same as the optimal semicoupling between A„ and fi n and 
similarly the optimal semicoupling between A and fx is the same as the optimal semicoupling between A 
and fi. Moreover, because for any n A„ is compactly supported with total mass N the vague convergence 
X n — > X implies weak convergence of X n — > A. 

iii) Now we are in a setting where we can apply the previous Lemma. Set K = 2R 2 and define p, n , p, as 
above. Then q n and q are optimal couplings between X n and fl n and A and fx respectively wrt to the cost 
function c(-, •). The cost function c is continuous and M and M U {3} are Polish spaces. Hence, we can 
apply the stability result of the classical optimal transportation theory (e.g. Theorem 5.20 in [Vil09] ) to 
conclude that q n — > q weakly and therefore q n — > q weakly. □ 

Take a pair of equivariant random measure (A*, A**) with A w <C m as usual. For a given w € 11 we want 
to apply the results of the previous Lemma to a fixed realization (A",/i w ). Then, for any bounded Borcl 
set A C M, there is a unique optimal semicoupling q M A between A" and 1am"j that is, a unique minimizer 
of the cost function Cost among all semicouplings of A" and 1aM w - 

Lemma 3.10. For each bounded Borel set A C M the map is measurable. 

Proof. We saw that the map $ : (A", 1aM w ) = 9a is continuous. By definition of random measures the 
map id h-> (A", IaA 4 ") is measurable. Hence, the map 

u^$(A",l Af x u ) = ct% 

is measurable. □ 
The uniqueness and measurably of q^ allows us to finally deduce 

Theorem 3.11. (i) For each bounded Borel set A C M there exists a unique semicoupling Qa 0/ A*P 
and (l^/x*)P which minimizes the mean cost functional £o5t(.). 

(ii) The measure Qa can be disintegrated as Q A(dx , dy , did) := q^dx, dy) P(dw) where for ¥ '-a. e. ui the 
measure q 1 ^ is the unique minimizer of the cost functional Cost(.) among the semicouplings of A" and 
1 \/r . 

(Hi) £ost(Q A ) = / a Cost(g^)P(da;). 

Proof. The existence of a minimizer is proven along the same lines as in the previous proposition: We 
choose an approximating sequences Q n in M.(M x M x f2) - instead of a sequence q n in A4(M x M) 
- minimizing the lower semicontinuous functional Cost(-). Existence of a limit follows as before from 
tightness of the set of all semicouplings Q with £ost(Q) < 2infg <£ost(Q). 
For each semicoupling Q of A* and IaM* with disintegration as q'F we obviously have 

<tosi(Q)= [ Cost( g w )dP(w). 
Jn 

Hence, Q is a minimizer of the functional £ost(-) (among all semicouplings of A* and IaM*) if an( i onr y 
if for P-a.e. u> £ f2 the measure q u is a minimizer of the functional Cost(.) (among all semicouplings of 
A" and 1am"0- 

Uniqueness of the minimizer of Cost(-) therefore implies uniqueness of the minimizer of Cost(-). □ 
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Corollary 3.12. The optimal semicouplings Qa = o'a^ a?~e equivariant in the sense that 

QgA(gC,gD,9 g u>) — Qa(C, D,uj), 

for any g G G and C,D G B(M). 

Proof. This is a consequence of the equivariance of A* and p* and the fact that q^° A is a deterministic 
function of A" and 1aH u '■ □ 

4. Uniqueness 

The aim of this section is to prove Theorem the uniqueness of optimal semicouplings. Moreover, the 
representation of optimal semicouplings, that we get as a byproduct of the uniqueness statement, allows 
to draw several conclusions about the geometry of the cells of the induced allocations. 
Throughout this section we fix two equivariant random measures A* and p* of unit resp. subunit intensi- 
ties on M with finite optimal mean transportation cost c etOC . Moreover, we assume that A* is absolutely 
continuous. 

Proposition 4.1. Given a semicoupling q w of A" and p w for fixed uj G then the following properties 
are equivalent. 

(i) For each bounded Borel set A C M, the measure ImxAQ u * s the unique optimal coupling of the 
measures A^(-) := q u (-,A) and IaH? ■ 

(ii) The support of q" is c- cyclically monotone, more precisely, 

N N 

^c{xi,yi) < ^c(x l ,y l+1 ) 

i=l i=l 

for any JVeN and any choice of points (xi, j/i), . . . , (xjv,2/at) in supp(g^) with the convention 

VN+l = J/1- 

(hi) There exists a nonnegative density p u and a c-cyclically monotone map T u : {p" > 0} — > M 
such that 

(4.1) g w = (Jd,T")Jp"A w ). 

Recall that, by definition, a map T is c-cyclically monotone iff the closure of its graph {(x, T(x)) : 
x G {p u > 0}} is a c-cyclically monotone set. 



Proof. {Hi) =>■ (it) =>■ (t) follows from Theorem 3.1 and Lemma 3.2 



(i) =>■ (tit) : Take a nested sequence of convex sets (K n ) n such that K n /*■ M . By assumption 1mxK„Q 
is the unique optimal coupling between A^- <C m and lfc n pL UJ . By Proposition 



3.6 



or Theorem 



3.1 



there exists a density p£ and a map T% : {p% > 0} -> M such that lMxK n <t = (^fn)*(Pn^K n )- Set 
— {p 1 ^ > 0}. As K n C K n+ i we have A£ C 4" +1 . Subtransports of optimal transports are optimal 
again. Therefore, we have 

implying p^+i = Pn on A^. Hence, the limits 

T u = limT" A u = limA£ and p u = lim^ 

n n n 

exist and define a c-cyclically monotone map T u : A u — > M such that on A u x M : 

□ 

Remark 4.2. In the sequel, any transport map T u : A u — > M as above will be extended to a map 
T w : M -> M U {3} by putting T""(x) := 3 for all x G M \ A" where 3 denotes an isolated point added 
to M ('point at infinity', 'cemetery'). Then ( |4.1[ ) reads 

(4.2) = (W, T") + p^A" on M x M. 

Moreover, we put c(x,T UJ (x)) — c(x,9) := for x G M \ A". If we know a priori that p u (x) G {0, 1} 
almost surely (4.2) simplifies to 

(4.3) <f = (Id, r), A" on M x M. 

Definition 4.3. A semicoupling Q = q'V of A* and p* is called locally optimal iff some (hence every) 
property of the previous proposition is satisfied for P-a.e. u> G £1. 
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Figure 3. The left picture is a semicoupling of Lcbcsgue and 36 points with cost func- 
tion c(x,y) — \x — y| 4 . In the right picture, the five points within the small cube can 
choose new partners from the mass that was transported to them in the left picture 
(corresponding to the measure A 1 ^). If the semicoupling on the left hand side is locally 
optimal, then the points in the small cube on the right hand side will choose from the 
gray region exactly the partners they have in the left picture. 



Remark 4.4. (i) Asymptotic optimality is not sufficient for uniqueness and it does not imply local 
optimality: Consider the Lebesgue measure A* = C = A and a Poisson point process \i* of unit intensity 



on W l . Let us fix the cost function c(x,y) = \x — y\ 2 . Lemma 2.14 shows that Co the optimal mean 
transportation cost on the unit cube, that is the cost of the optimal semicoupling of A and lm i)<*/U*, is 
strictly less than the optimal mean transportation cost on a big cube, say [0, 10 10 ) d . Moreover, for any 
semicoupling q' between A and p' Lemma 2.14 implies that 

£(<?•) = liminf — — E[Cost(l R< i X B w O] , 

n— i-oo A(-D n ) 

for B n = [— 2 n_1 , 2™~ 1 ) d . In other words, it is more costly to transport in one big cube than in many 
small cubes separately. For all n <G N, let p n : ft x R d — > [0, 1] be the unique optimal density for the 
transport problem between A and ls„/i*, that is the optimal semicoupling is given by g* = (id, T*)*(p* A). 
Let k* be the following semicoupling between A and 

K^idx^dy) = q^(dx,dy) + a^(dx,dy), 

where er^ is the unique optimal coupling between (p% — p^) ■ A and ls n \B /x w . Let /„ : n x R d ->■ [0, 1] be 
such that l_Rd x (B„\-B Q )'J'n — (i<l,^) t (/^). Denote by W2 the expectation of the usual L 2 — Wasserstein 
distance. Then, we can estimate using the triangle inequality 

£ost 1/2 «) = €ost 1/2 (q')+W 2 ((pn- Po) ■X,l {Bn \B )^') 

< £ost 1/2 (q' ) + €ost 1/2 (q' n ) + W 2 ((p„ - pa) ■ XJn ■ A). 

Set Zi — p,*(Bi). Note that (p„ — Po) and coincide on a set of Lebesgue measure of mass at least 
Z% — Zq. This allows to estimate W2((/) m — po) • A, f n ■ A) very roughly from above (for similar less rough 
and much more detailed estimates we refer to the cost estimates in section 5 of [HS10 ). We have to 
transport mass of amount at most Zq at most a distance R\ = 2h ■ z}J d + 2\fd2 n for some constant h, e.g. 
h = (r(| + l)) 1 ^ would do. Indeed, must be supported in a h ■ Z x l d neighbourhood of B n because 
we could otherwise produce a cheaper semicoupling (see Lemma 3.9). This gives using the estimates on 
Poisson moments of Lemma 5.11 in [HSlOj 

Wl((p n -p )-X,f n -X)<E[Rl-Z ] 

< C\ (E[Z 2 ]^ 2 E[Z^ d }^ 2 + X{B n ) 2 ' d + 2 • E[Z 2 ] 1 / 2 A(S„) 1 / d E[Z2/ rf ] 1 /2 

< C 2 X{B n ) 2 ' d , 

for some constants C\ and C 2 . In particular, if we take d > 3 this shows that 

liminf T7~R~^^ osi ( K '^ = liminf T7^ £osi ^*)- 
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Hence, as in Proposition 2.20 we can show that converges along a subsequence to some semicoupling 
k' between A and /i* which is asymptotically optimal but not locally optimal. 

(ii) Local optimality does not imply asymptotic optimality and it is not sufficient for uniqueness: For 
instance in the case M = M. d , c(x, y) = \x — y\ 2 , given any coupling q* of C and a Poisson point process 
ff and z € R d \ {0} then 

F{dx,dy) :=q u] (d(x + z),dy) 
defines another locally optimal coupling of C and /i* . Not all of them can be asymptotically optimal. 

(iii) The name local optimality might be misleading in the context of semicouplings. Consider a Poisson 
process fi* of intensity 1/2 and let q' be an optimal coupling between l/2£ and /iV Then, it is locally 
optimal (see Theorem 4.6) according to this definition. However, as we left half of the Lebesgue measure 
laying around we can everywhere locally produce a coupling with less cost. In short, the optimality does 
not refer to the choice of density only to the use of the chosen density. 

(iv) Note that local optimality — in contrast to asymptotic optimality and equivariance — is not 
preserved under convex combinations. It is an open question if local optimality and asymptotic optimality 
imply uniqueness. 

Given two random measures 7*, rf : O — > M.{M) with 7 W (M) = 77" (M) < 00 for all uj € we define the 
optimal mean transportation cost by 

€ost(Y,v') ■= inf {£ost(q*) : q u € H{-f ,rf) for a.c. wefi}. 

Given a semicoupling q* of A* and /x* and a bounded Borcl set A c M recall the definition of A^ from 
Prop. 4.1 We define the efficiency of the semicoupling q* on the set A by 

€osi(\ A , IMS) 

e ^ } := ^aw) ' 

It is a number in (0, 1]. The semicoupling q* is said to be efficient on A iff zff A (q*) = 1. Otherwise, it 
is inefficient on A. As noted in the remark above in the case of true semicouplings this notion might 
mislead the intuition. 

Lemma 4.5. (i) q* is locally optimal if and only if cff A (q*) = 1 for all bounded Borel sets A C M. 
(ii) eff j4 (g*) = 1 for some Ad M implies tff A i{q') = 1 for all A' C A, where we set 0/0 = 1. 

Proof, (i) Let A be given and uj £ Vt be fixed. Then ImxAQ 1 * 1 is the optimal semicoupling of the measures 
X A and 1aH u if and only if 

(4.4) Cost(l M xA<n = Cost {X% , l A fj.") . 

On the other hand, cff A (q*) — 1 is equivalent to 

E [Cost(lMxA<Z*)] = E [Cost (X' A , l Al j, m )] . 



The latter, in turn, is equivalent to (4.4) for P-a.e. ui £ 1L 

(ii) If the transport q* restricted to M x A is optimal then also each of its sub-transports. □ 

Remember that due to the stationarity of P equivariance of q* translates into invariance of its distribu- 
tion. The next Theorem is a key step in establishing uniqueness because it shows that every optimal 
semicoupling is induced by a map. 

Theorem 4.6. Every optimal semicoupling between A* and ji* is locally optimal. 

Proof. Assume we are given a semicoupling q' of A* and fi' that is equivariant but not locally optimal. 
According to the previous lemma, the latter implies that there is g € G and r £ N such that q* is not 
efficient on gB r7 i.e. 

V = ^ 9 bM') < L 

By invariance, this implies that 77 = eff hB (<?*) < 1 for all h € G. Hence, for any h S G there is a coupling 
q' B of A* B and lhB r fJ-*j the unique optimal coupling, which is more efficient than lMxhB r q* , i-e. such 
that 

E[Cost(<7* Br )] < n-E[Cost(l MxhBr q')}. 
Moreover, because of the equivariance of q* we have q^ B (dx, dy) — qJ^B (d(9 x )i d(gy)) (see also Corollary 



3.12). Hence, all convex combinations of the measures qf lB will have similar equivariance properties. 
We would like to have the estimate above also for the restriction of q' iB to M x hB in order to produce 
a semicoupling with less transportation cost than q*. This is not directly possible as we cannot control 
the contribution of hBo to the cost of <j* B . However, we can use a trick which we will also use for the 
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construction in the next section that will give us the desired result. Remember that A r denotes the 2 r 
neighbourhood of the identity in the Cayley graph of G. Set 

QjiBo = TT i ^My.hB Q q*gB r - 
|Arl g£hA r 



Then we have 



£°st(q' hBo ) = r V E [ c(x,y)q' gBr (dx,dy) 

l A H „ c u UMxhBo 



g£hA 

|A r | 



MxhB 

c ( x ,y)VhB r {dx,dy) 



MxhB r 



= r^T-E[C0St(^ flr )]. 

The second equality holds because fixing hB$ and summing over all gB r containing hBo is, due to the 
invariance of q* Br , the same as fixing hB r and summing over all gB$ contained in hB r . Put q* — 
^2heG 9*B ' ^ cons t ruc ti° n ; Q* is an equivariant semicoupling of A* and /i*. Furthermore, for any g G G 
we have 

E[CoSt(<7^g )] < ?7 • E[lMxgB 9*]- 

This means, that q* is not asymptotically optimal. 

□ 

Remark 4.7. We really need to consider r > 1 in the above proof as it can happen that q* is efficient 
on every fundamental region but not locally optimal. Indeed, consider p, = J2 Z £i, 2 and let q' denote 
the coupling transporting one quarter of the Lebesgue measure of the square of edge length 2 centered at 
z to z. This is efficient on every fundamental region, which contains exactly one z £ Z 2 , but not efficient 
on say [0, 5) 2 . 

Theorem 4.8. Assume that p,' has intensity one, then there is a unique optimal coupling between A* 
and p* . 

Proof. Assume we are given two optimal couplings q' and q°- Then also q' :— ^q' + ^q' is an optimal 
coupling because asymptotic optimality and equivariance are stable under convex combination. Hence, 
by the previous theorem all three couplings - q*, q* and q* - are locally optimal. Thus, for a.e. to by the 
results of Proposition |4. I| there exist maps T" ,T£ ,T U such that 

$T"(x){dy) A" (eta) = q u (dx,dy) 

~5 Tr{x) (dy) + ~S T u, {x) (dy) J \"(dx) 



This, however, implies T"(x) = T£(x) for a.e. x € M. Thus qf — q% ■ (By Lemma 2.22 we know that 
every invariant semicoupling between A* and \i* has to be a coupling.) □ 



Before we can prove the uniqueness of optimal semicouplings we have to translate Proposition 3.5 to this 
setting. 

Proposition 4.9. Assume p* has intensity /3 < 1 and let q* be an optimal semicoupling between A* and 
p* . Let (7Ti)*g* = p ■ A* for some density p : fi x M [0,1]. Then, 

<f = q* + (id x -p) ■ A*) 

is the unique optimal coupling between A* and p,' := p,' + (1 — p) ■ A*. 

Proof. Because q' is equivariant by assumption also pA*(-) = q*(-,M) is equivariant. But then p,' — 
p,' + (1 — p) ■ A* is equivariant. Moreover, by assumption we have (£((?") = £(<?*) < 00 which implies 

inf £(«*) < 00. 
K*en e (A , ,A*) 

By the previous theorem, there is a unique optimal coupling re* between A* and /t* given by k* — 
(id,S)*\*. Moreover, 

£(«•) < £(«•) = £(«•)• 

Because 5* A* = //* there is a density / such that • A*) = (1 — p) ■ A*. Indeed, for any g £ G we can 
disintegrate 

lMx 9Bo ^(dMy) - ^• 9 (^)(M w (<iy) + (l-p a, (y))A w (dy)). 
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The measure J2 g <£G K y' 9 (d x )((l ~ P u {y))^ u {dy)) does the job. In particular this implies that 

ft* = (id x S%((1 - /) • A*) 

is a semicoupling between A* and p*. The mean transportation cost of ft* are bounded above by the 
mean transportation cost of k* as we just transport less mass. Hence, we have 

£(ft') < £(*•) < <t(q') = 

As q* was assumed to be optimal, hence asymptotically optimal, we must have equality everywhere. By 
uniqueness of optimal couplings this implies that q* = k' almost surely. □ 

Lemma 4.10. Assume p* has intensity @ < 1 and let q* = (id, T%(p ■ A*) be an optimal semicoupling 
between A* and p* . Then, on the set {0 < p^ < 1} we have T u (x) = x. 

Proof. Just as in the previous proposition consider q* — (id, if?)*A* the optimal coupling between A* and 
fx*, q* is concentrated on the graph of S and therefore also q* has to be concentrated on the graph of S. 
In particular, this shows that S = T almost everywhere almost surely (we can safely extend T by S on 
{p = 0}). But on {p < 1} we have S(x) = x. Hence, we also have T(x) = x on {0 < p < 1}. □ 

This finally enables us to prove uniqueness of optimal semicouplings. 

Theorem 4.11. There exists at most one optimal semicoupling of A* and p* . 

Proof. Assume we are given two optimal semicouplings q\ and q°. Then also q' := \q° + is an 
optimal semicoupling. Hence, by Theorem 4.6 all three couplings - q', q' and q* - are locally optimal. 
Thus, for a.e. uj by the results of Proposition 4.1 there exist maps Tf ,T£ , T w and densities , P2 , 
such that 

6 T *( x) (dy)p"(x)\ u (dx) = q lJ (dx,dy) 

\sT nx) {dy)pt(x) + \8 Tnx) (dy)pt(x) \ X»(dx) 

This, however, implies Tf(x) = T%{x) for a.e. x G {p^ > 0} n {p^ > 0}. In particular, all optimal 
semicouplings are concentrated on the same graph. To show uniqueness, we have to show that p^ = 
almost everywhere almost surely. To this end, put = {p^ > p%}. Assume X bJ (A^) > 0. On A" we 
have p u < 1. Hence, by the previous Lemma we have T w (x) = T"(x) — Ttffa) — x. Similarly, on 
^2 = {/°2 > Pi} we nave T u (x) = x. Hence, we have 

(T u )- l (A%)f\A% = 0, 

because Ai C (T u )^ 1 (Ai) for i = 1,2. As q^ and q% are semicouplings, we must have p u (A) = pf ■ 
A w ((T ,u )~ 1 (vl)) for i = 1, 2 and any Borel set A. Putting this together gives 

H U (A%) = pff ■ A"((T")- 1 (^)) 

= pi ■ X U ((T U )- 1 (A^) n A") + ft ■ x u ((T u y 1 (A u 1 ) n {p^ = p^}) 
> P % ■ x^T")- 1 ^) n A") + ft ■ A-CCT")- 1 ^) n K = P2 }) 

This is a contradiction and therefore proving X UJ (Ai) = 0. Thus, pf = p% almost everywhere almost 
surely and q' = q*. □ 

4.1. Geometry of tessellations induced by fair allocations. The fact that any optimal semicoupling 
is locally optimal allows us to say something about the geometries of the cells of fair allocations to point 
processes. The following result was already shown for probability measures in section 4 of [Stullaj and 



also in |AHA92| . We will use the representation of optimal transportation maps recalled in Remark 3.4 



Corollary 4.12. In the case i?(r) = r 2 , given an optimal coupling q* of Lebesgue measure C and a point 
process p' of unit intensity in M = M. d (for a Poisson point process this implies d > 3 as otherwise the 
mean transportation cost will be infinite, see Theorem 1.3 in [HSlOj ) then for a.e. ui G f2 there exists a 
convex function ip u : R d —> R ( unique up to additive constants ) such that 

In particular, a 'fair allocation rule' is given by the monotone map T u — V<^ w . 
Moreover, for a.e. uj and any center £ G S(cj) := supp(/Lt w ), the associated cell 
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is a convex polytope of volume /x w (C) G N. If the point process is simple then all these cells have volume 
1. 

Proof. See Corollary 3.10 of [HSIO] . □ 

Corollary 4.13. In the case d(r) = r, given an optimal coupling q* of m and a point process /i* of unit 
intensity on M with dim(M) > 2, there exists an allocation rule T such that the optimal coupling is given 
by 

Moreover, for a.e. uj and any center £ G 2(cj) := supp(/x^), the associated cell 
is starlike with respect to C- 



Proof. By Proposition 4.1 we know that = limn._j.oo T£ , where T£ is an optimal transportation map 
from some set to K n . From the classical theory (see |Bre91( IGM961 ) we know that, 

T%(x)=£o & -d(x,^ ) + h a > -d(x,0 + h € H(w) n K n , ^ Co- 
Hence, the cell can again be written as the intersection of "halfspaces" := {x : —d(x,(_o) + > 
—d(x,^j) + &£.}. Therefore, it is sufficient to show that for any z G H° the whole geodesic from z to Co 
lies inside For convenience we write &o(x) = — Co) + ^ an d ^jO^) = —d{x,^j) + b^. 
Assume Co G d-ff? and w.l.o.g. 6^ = 0. Then, we have 

*o(eo) = o = $ i (e )=^6c = d(( 3 -,&). 

The set N — {z G M : d(^j,z) = d(£j,£o) + d(£o, z)} is a m-null set. For all z ^ _V we have 

$j(z) = -d(£ jt z) + b^ > -d(^,^o)+b Cj - d(£ ,z) = $„(*) 

This implies that m(T~ 1 (£j)) = contradicting the assumption of T being an allocation. Thus, Co 4- 
and in particular T(£ ) G S = supp(/i). 

Assume T(£ ) ^ Co- Then, there is a ^ ^ Co such that T(£ ) = G, i.e. $j(Co) = -d(£o,&) + b^ > b 5o = 
$o(Co)- Then, we have for any p G M,p ^ Co 

-d(p, Cj) + h 3 > - d (P, Co) - d(Co, Cj) + b« 3 > ~d(p, Co) + b £o . 
This implies, that m(T _1 (Co)) = contradicting the assumption of T being an allocation. Thus, T(£o) = 
T„(Co)=Co. 

Take any w G T n x (Co) (hence, <I > o(w) > $j(«.) for all j ^ 0) and p G M such that d(£o,w) = d(£ ,p) + 
d(p,w), i.e. p lies on the minimizing geodesic from Co to w. Then, we have for any j ^ by using the 
triangle inequality once more 

-d(p>£o) + &£ = -d(£o,to) + d(p,iu) + b io 

> -^(Co, w) + b ?0 + d(w, Cj) - d(p, Cj) 

> -rf(p,Cj)+b^, 

which means that $o(p) > ^j(p) f° r au j 7^ 0. Hence, p G -H^ proving the claim. □ 

Remark 4.14. i) Questions on the geometry of the cells of fair allocations are highly connected to the 
very difficult problem of the regularity of optimal transportation maps (see [MTW051 ILoe09[ IKM10] ) . 
The link is of course the cyclical monotonicity. The geometry of the cells of the "optimal fair allocation" 
is dictated by the cyclical monotonicity and the optimal choice of cyclical monotone map to get an 
asymptotic optimal coupling. 

Consider the classical transport problem between two probability measures one being absolutely contin- 
uous to the volume measure on M with full support on a convex set and the other one being a convex 
combination of N Dirac masses. Assume that the cell being transported to one of the N points is not 
connected. Then, it is not difficult to imagine that it is possible to smear out the Dirac masses slightly 
to get two absolutely continuous probability measures (even with very nice densities) but a discontinuous 
transportation map. 

ii) Considering IP cost on M. d with p {1,2}, the cell structure is much more irregular than in the two 
cases considered above. The cells do not even have to be connected. Indeed, just as in the proof of the 
two Corollaries above it holds also for general p that T u (x) = Co iff &o{x) > $i(x) for all i =t where 
$i(x) — — \x — C-I p + b. for some constants bj (see also Example 1.6 in [GM5B]). By considering the sets 
<&i = $0 it is not difficult to cook up examples of probability measures such that the cells do not have to 
be connected. 
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In the case that p € (0, 1) similar to the case that p = 1 we always have that the center of the cells lies in 
the cell, that is T(£j) = & for all € supp(^*) because the cost function defines a metric (see |GM96j ). 

hi) As was shown by Loeper in section 8.1 of [Loc09 the cells induced by the optimal transportation 
problem in the hyperbolic space between an absolutely continuous measure and a discrete measure with 
respect to the cost function c(x,y) — d 2 (x,y) do not have to be connected. In the same article he shows 
that for the same problem on the sphere the cells have to be connected. In [vN09 von Nessi studies more 
general cost functions on the sphere, including the LP cost function c(x,y) = d p (x,y). He shows that in 
general for p ^ 2 the cells do not have to connected. This suggests that on a general manifold the cell 
structure will probably be rather irregular. 

5. Construction 

We fix again a pair of equivariant random measures A* and /x* of unit resp. subunit intensity with finite 
optimal mean transportation cost c ei00 such that A* is absolutely continuous. Additionally, we assume 
that G satisfies some growth condition. Recall that the 2 r neighbourhood of the identity element in the 
Cayley graph A(G, S) of G is denoted by A r and the range of its action on the fundamental region by 
B r , i.e. B r = U g eA S^o- Then, we assume that for any g G G 

(5.1 1 , t , 1 -> as r -> oo, 

|A r | 

where \A\ denotes the cardinality of A. This of course implies for any g € G 

m(B r AqB r ) 

y , ' -> as r -> oo. 
m(B r ) 

The aim of this section is to construct the optimal semicoupling and thereby proving Theorem |1.2| The 
construction is based on approximation by semicouplings on bounded sets. We will also show a nice 



convergence result of these approximations, proving Theorem 1.3 The proofs in the first part of this 
section can mostly be copied from the respective results in [HS10 . Therefore, we omit some of them and 
only stress those where something new happens. 

5.1. Symmetrization and Annealed Limits. The crucial step in our construction of optimal semi- 
couplings between A* and fj,' is the introduction of a symmetrization or second randomization. We want 
to construct the optimal semicoupling by approximation of optimal semicouplings on bounded sets. The 
difficulty in this approximation lies in the estimation of the contribution of the fundamental regions gB 
to the transportation cost, i.e. what does it cost to transport mass into gBo? How can the cost be 
bounded in order to be able to conclude that the limiting measure still transports the right amount of 
mass into gB$l The solution is to mix several optimal semicouplings and thereby get a symmetry which 



will be very useful (see proof of Lemma 5.1 (i)). One can also think of the mixing as an expectation of 
the random choice of increasing sequences of sets hB r exhausting M. 

For each g £ G and r G N, recall that Q g B r denotes the minimizer of Cost among the semicouplings 



of A* and l g B r IJ'' as constructed in Theorem 2.13 It inherits the equivariance from A* and /i*, namely 

QgA(g-,g-,0c 

QhgB r - Put 



Q g A{g', g',0gLu) = Qa('> (see Corollary 3.12). In particular, the stationarity of P implies (Th)*Q g Br = 



1 

I 

h£gA r 



Q r g (dx,dy,duj) := l g B (y)r^j ^2 QhB r ( dx > d y> duj )- 



The measure Q g defines a semicoupling between A* and l g s /z*. It is a deterministic, fractional allocation 
in the following sense: 

• for any uj it is a deterministic function of A w and /i w and does not depend on any additional 
randomness, 

• for any to the first marginal is absolutely continuous with respect to A" with density < 1. 

The last fact implies that the semicoupling Q r g is not optimal in general, e.g. if one transports the 
Lebesgue measure to a point process. The first fact implies that all the objects derived from Q r g in the 
sequel - like and Q°° - are also deterministic. Moreover, Q r g shares the equivariance properties of 
the measures QhB r - 

Lemma 5.1. (i) For each r E N and g e G 

c(x 1 y)Q r (dx,dy 1 duj) < c^. 



Mx g B xU 
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h 3 B r 






h\B r 





h,2B r 

FIGURE 4. Schematic picture of the mixing procedure. 



(ii) The family (Q g ) r eN of probability measures on M X M X f2 is relatively compact in the weak 
topology. 

(iii) There exist probability measures and a subsequence (r;); S N such that for all g £ G: 

Q r „ l — > weakly as I — > oo. 

Proof, (i) Let us fix g £ G and start with the important observation: For given r £ N and g £ G 
averaging over all hA r with h £ gA r has the effect that "g attains each possible position inside A r with 
equal probability" (see also the proof of Theorem 4.6 1. 
Hence, together with the invariance of QkB r we obtain 



c(x, y)Q r g {dx, dy, dw) 

MxgB xU 

r-f V" / c(x,y)Q hBr ( dx idy,du) 

V H L _ . JMxaBnxn 



MxgB r xU 



1 MxgB xQ. 

c(x,y)Q gBr (dx,dy,du]) 



|A., , 

1 1 hegA 

1 

= \h C9Br 

by definition of Coo. 

(ii) In order to prove tightness of (QpreNi let (gBo)i denote the closed 1-neighborhood of gB in M. 
Then, 

Q r g(E(9 B o)l,gB ,£l) < / c{x,y)Q r g (dx,dy,dui) 

UK 1 ) JMxgBoxQ. 
1 

- m Coc ' 

Since i?(i) — > oo as 2 — > oo this proves tightness of the family (Qp r gN on M x M x f2. (Recall that f2 
was assumed to be compact from the very beginning.) 

(iii) Tightness yields the existence of and of a converging subsequence for each g £ G. A standard 
argument ('diagonal sequence') then gives convergence for all g £ G along a common subsequence (G is 
countable as it is finitely generated). □ 

Note that the measures inherit as weak limits the property Q^ g {h-,h-,6h-) — Q^°(-,',-) from the 
measures Q r g (see also the proof of the equivariance property in Proposition 2.20). The next Lemma 
allows to control the difference in the first marginals of and for g ^ h. This is the first point 
where we use the growth condition. 
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Lemma 5.2. i) For all I > there exists numbers e r (l) with e r (l) — > as r — > oo s.t. for all 

g, g € G and all r N 

QhbM)<tt1 QhB r (A) + e r (d A {g,g'))- sup Q hBr (A) 

|Ar| hes'A,. |VV| hegK h&g' A r 

for any Borel set A C M x M x O. 
ii) For all gi, . . . , g n eG, ail r e N and all Borel sets A C M, D a Q 

n / n \ 

^Q;,(AM,fl) < l + ^e r (d A ( 5l , 9 i)) ■ A(-D,A), 

i=l \ i=l / 

where \{D,A) := J D J A X u (dx)f(du). 
Proof, (i) First note that for all g, g € G and rdNwe have 

g e gA r g e g A r . 

In this case, for h S gA r with g € ftA r we also have h £ g A r and g € ftA r . Moreover, 

€ gA r : g $ hA r }\ , 

^ ^a&s )), 

for some e r (l) with e r (Z) — > as r — > oo. One possible choice for e r is 

/j (i , » |A r A. g A r | 
e r (d A {id,g)) = — , 

which tends to zero as r tends to infinity for any g ^ id by assumption. This implies that for each pair 
g,g G G and each r G N 

e gA r : </ e /iA r }| 



> 1 - e r {d A {g,g )). 



|A r | 

Therefore, for each Borel set A C M x M x!! 

ro QhB r {A)<^- V g, iSr (A)+e r (d A (5,5'))- sup Q, iBr (A). 

|Ar| hes 'A r 1 r| ^e 9 A. heg'A r 

(ii) According to the previous part (i), for each Borel sets A C M, D C il 

n 

i=i 

n 

= Eul E Q hBr (A9iB ,D) 

i=l I r l h£ gi A r 



n I 

E tat E QhB^gi^D) 

i=l V I r l hEgiA r 

< (l + f2^(dA(9i,9i)))MD,A) 



- /Am /] QhB r (A,giB ,D) + e r (dA{gi,gi))- sup Q hBr .( A i 9i B o, D) 

hegiA,, 



□ 



Having these results at our hands we can copy basically line to line the respective proofs from HS10 
(Theorem 4.3 and Corollary 4.4) to get 

Theorem 5.3. The measure Q°° := X^eG * s an optional semicoupling of A* and /z*. 

Corollary 5.4. (%) For r — > oo, i/ie sequence of measures Q r :— ^ Q^, r G N, converges vaguely to 

i/ie unique optimal semicoupling Q°° . 
(ii) For each g S G and r € N pitf 

Q r g {dx,dy,duj) := -j^— ^ QhB r { dx , d V, duj )- 
The sequence (Q 7 Q ) r eN converges vaguely to the unique optimal semicoupling Q°° . 
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Corollary 5.5. Denote the set of all semicouplings of A* and fi* by II 

1 



Then it holds 



inf lim inf 
g'en 3 r->oo m(B r ) 

= lim inf inf 



E 



c(x,y)q*{dx,dy) 



1 



In particular, we have 



9'en, m(B r ) 



inf £(g*) 
9*en, 



MxB. 



E 



Afxj 



c(x,y)q'(dx,dy) 



inf C(g') 

<?*en es 



Proof. For any semicoupling q* we have due to the supremum in the definition of <£(•) that 

1 



lim inf 



m(B r 



-E 



c(x,y)q'{dx,dy) 

IMxB r 

Hence, the left hand side is bounded from above by inf g . e n s <£(?*)■ However, we just constructed an 
equivariant semicoupling, the unique optimal semicoupling Q°° which attains equality, i.e. with Q°° = q'V 

1 



lim inf ■ 



-E 



c(x,y)q*(dx,dy) 



MxB r 



= C(Q°°). 



m(B r ) 

Hence, the left hand side equals inf 9 . e n s £(<?*)• 
The right hand side equals lim inf r-^oo c r which is bounded by c c 
our construction, the asymptotic transportation cost of Q°° are bounded by the right hand side, i.e. 

£(Q°°) < lim inf c r 



inf g . 6 n s by Lemma 



2.14 



by Lemma 5.1 Hence, also the right hand side equals inf g . e n s £■(?*)■ Thus, we have equality. 



By 



□ 



Remark 5.6. i) Because of the uniqueness of the optimal semicoupling the limit of the sequence Q r does 
not depend on the choice of fundamental region. The approximating sequence (Q r ) r eK does of course 
depend on B and also the choice of generating set S that defines the Cayley graph. 

ii) In the construction of the semicoupling Q°° we only used finite transportation cost, invariance of Qa 

in sense that (Th)*QA = QhA an d the amenability assumption on G. The only specific property of A* and 
/i* that we used is the uniqueness of the semicoupling on bounded sets which makes is easy to choose a 
good optimal semicoupling Q g B r ■ Hence, we can use the same algorithm to construct an optimal coupling 
between two arbitrary random measures. In particular this shows, that = c e ^ (see also Proposition 



2.20). 



Indeed, given two arbitrary equivariant measures v* and fi° of unit respectively subunit intensity. For 
any r £ N let Qb t — q% V he an optimal semicoupling between v* and IaM*- In particular, we made 
some measurable choice of optimal semicoupling for each u) (they do not have to be unique), e.g. like in 
Corollary 5.22 of |Vil09j . Define Q g B r via q e g 3 g r (d(gx),d(gy)) := qg r (dx,dy). Due to equivariance, this is 

again a measurable choice of optimal semicouplings. Stationarity of P implies (t^)*Qb,. = Qkb t - Hence, 
by the same construction there is some optimal semicoupling Q°° of v* and fi* with cost bounded by Cqq. 

5.2. Quenched Limits. According to section |1J the unique optimal semicoupling between A* and fi' 
can be represented on M x M x Q as 

Q°°(dx,dy,duj) = S Tix<u) (dy) p"(x)\"(dx) P(du) 

by means of a measurable map 

T : M xil^ MU{9}, 

defined uniquely almost everywhere and a density p u . Similarly, for each g £ G and rgN there exists a 
measurable map 

T gtr : M x fi ->• M U {9} 
and a density such that the measure 

Q gBr .{dx,dy,(kS) = S Tg Ax ^(dy) Pg r \ u (dx)¥(duj) 
on M X M X fi is the unique optimal semicoupling of A* and l g B r A*" • 
Proposition 5.7. For every g 6 G 

Tg tT (x,ui) — > T(x,uS) as r — > oo in A* ®P -measure. 

The claim relies on the following two Lemmas. For the first one we use the growth assumption once more. 
The second one is a slight modification (and extension) of a result in [A"mb03 . 
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Lemma 5.8. i) Fix uj 6 Q and take two disjoint bounded Borel sets A, B C M. Let = 

(id,T^)^(p L ^X UJ ) be the optimal semicoupling between A" and 1aP- w - Similarly, let q B and q%jg 
be the unique optimal semicouplings between A" and \bP? respectively \a\jbP^ with transport 
maps Tg and T^ uB and densities p B and p^uB- Then, it holds that 

Paub( x ) > ^^{Pa( x )>Pb{ x )} X u a.s.. 

ii) For any g £ G and r £ N we have Pg tr (x) < p u (x) (A* ® P) a.s.. 

iii) For any g £ G we have lim,._ i . 00 p^ T (x) /* p u (x) (A* ® P) a.s.. 

Proof, i) Firstly, note that if {p^ > 0}Pi{p B > 0} = we have p%j B — Pa + Pb- Because of the symmetry 
in A and B it is sufficient to prove that Paub — P%- The proof is rather technical and involves an iterative 
choice of possibly different densities. 

For simplicity of notation we will suppress uj and write f = pb and h — paub and T — Tb, S — Taub- We 
will show the claim by contradiction. Assume there is a set D of positive A measure such that f(x) > h(x) 
on D. Put /_|_ := (/ — h) + and pi := T*(/_|_A). Let hi < h be such that S*(hiX) = pi, that is hi is a 
subdensity of h such that T*(/+A) = S^/iiA) (for finding this density we can use disintegration as in the 
proof of Proposition |4.9[ ) . 

If f{/i 1>0 }/i > / on some set Di of positive A measure, we are done. Indeed, as / is the unique Cost 
minimizing choice for the semicoupling between A and IbP the transport S*(lD 1 hiX) =: pi must be more 
expensive than the respective transport T t (ljj /+A) = pi for some suitable set D\. Hence, qAuB cannot 
be minimizing and therefore not optimal, a contradiction. 

If 1 {tnyojh < / we can assume wlog that T*(hiX) = pi and pi arc singular to each other. Indeed, if they 
are not singular we can choose a different hi because ls/i has to get its mass from somewhere. To be 
more precise, if h < hi is such that T*(h\) < pi we have T*((/+ + h)\) > pi. Therefore, there must be 
some density h! such that h! + h\ < h and S*((h' + hi)X) = + h)X). Because, /+ > on some set 

of positive measure and T*(/A) < S*(hX), there must be such an hi as claimed. 

Take a density h% < h such that S*(/i2^) = M2- If l{h 2 >o}^ > / on some set D2 of positive A measure, we 
are done. Indeed, the optimality of qs implies that the choice of / + and hi is cheaper than the choice of 
hi and h2 for the transport into pi + P2 (or maybe subdensities of these). 

If l{/i9>o}^- ^ / and {^2 > 0} H {/+ > 0} has positive A measure, we get a contradiction of optimality of 
qAuB by cyclical monotonicity. Otherwise, we can again assume that T^/^A) =: ps and P2 are singular 
to each other. Hence, we can take a density /13 < h such that S*(h,3\) = P3- 

Proceeding in this manner, because /+A(M) = hiX(M) > for all i and the finiteness of qs{M, M) one 
of the following two alternatives must happen 

• there is j such that l{h >o}/i > / on some set of positive A measure. 

• there are j 7^ i such that {hj > 0} n {hi > 0} on some set of positive A measure with / + = h . 

Both cases lead to a contradiction by using the optimality of qs, either by producing a cheaper semicou- 
pling (in the first case) or by arguing via cyclical monotonicity (in the second case). 

ii) Fix ui,g and r. Denote the density of the first marginal of Q 1 ^ by £/;• It is a convex combination of 
p% 1 with h £ fAi. For h £ G with d(g, h) < n we have gA r C hA r+n . Hence, we have Pg r < p£ r+n by 
the first part of the Lemma. Therefore, the contribution of Pg r (x) to (g r+n (x) is at least the number of 
h £ G such that d(g, h) < n divided by |A r+ „|. Hence, 

1X1 pl r {x)<Q r+n {x). 



By the assumption (5.1 1 we have 



,. \KA r AA r \ 

llm TaH = °< 

1 — s-00 |A r | 



for any finite K C G. If we take K = {h : d(h, id) — r} we can conclude 

|A„ +r | |^A n AA„| 

■ . 1 < 1 + ' nrn 1 -> 1 as n -t 00. 

|A„| " |A„| 

Fix e > 0. If Pg . r > e + p u on some positive (A* ® P)- set, we have that Cg, r +ni x ) > P u ( x ) + e / 2 on 
some positive (A* <£> P)— set for all n such that > 1 — e/2, because Pg r < 1 and thus p u < 1 — e. 

Denote this set by A, so A C M x ft. Then, we have Q r g +n (A x M) > Q°°(Ax M) + e/2 for all n big- 
enough. However, this is a contradiction to the vague convergence of Q r g to Q°° which was shown in the 
last section. 
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iii) The last part allows to interpret p^ r as a density of (p"A w ) instead of as a density of A w . We will 
adopt this point of view and show that r converges to 1 (A* <E> P) a.s.. 

Assume that p^ r (x) < 7 < 1 for all r £ N. Moreover, assume that there is k £ G and seN such that 
Pk s( x ) > 7- Then there is a t £ N such that gA t D fcA s . The first part of the Lemma then implies that 
Pg t( x ) ^ Pk s( x ) > 7 which contradicts the assumption of p 1 ^ r (x) < 7. Hence, if we have p^ r (x) < 7 < 1 
for all r € N on a set of positive (A* ® P) measure we must have p% s {x) < 7 for all A; € G ans seN on 
this set. Denote this set again by A, A C Af x f2. As is a convex combination of the densities pf t r it 
must also be bounded away from 1 by 7 on the set A. However, this is again a contradiction to the vague 
convergence of Q r g to Q°° . 

□ 

Lemma 5.9. Let X, Y be locally compact separable spaces, 9 a Radon measure on X and p a metric on 
Y compatible with the topology. 

(i) For all n £ N let T n ,T : X — > Y be Borel measurable maps. Put Q n (dx, dy) :— Sx n ( x )(dy)9(dx) and 
Q(dx,dy) :— ST( x ){dy)0(dx). Then, 

T n — > T locally in measure on X Q n — > Q vaguely in A4(X x Y). 

(ii) More generally, let T and Q be as before whereas 



Q n (dx,dy):= / 5 Tn ( x ^ x , } (dy) 0'(dx') 9{dx) 
Jx> 

for some probability space (X',21',0') and suitable measurable maps T n : X x X' — > Y. Then 

Qn — > Q vaguely in Ai{X x Y) T n {x,x') — > T{x) locally in measure on X x X . 

For a proof we refer to section 4 of [HS10J 

Proof of the Proposition. Firstly, we will show that the Proposition holds for 'sufficiently many' g £ G. 
We want to apply the previous Lemma. Recall that 

Q r g -> Q°° vaguely on M x M x CI, 

where 

Q°°(dx,dy,du) = 6 TiXiU) (dy)p u (x)\ u (dx)F(du}) 

and 

Q r g {dx,dy,du}) = QhB r { dx ^ dy, dw) = S Th _ r{x) {dy)pf l r (x)X UJ (dx)F(duj), 

' r ' hegB,. ' r ' hegB r 

with transport maps T, T/j jr : M x — > M U {3} and densities p, ph, r '■ M xO — > M + . The Lemma above 
allows to interpret ph, r as density of the measure pA*. Fix k £ G and let 6^ be the uniform measure on 
kA r . Take = p\° ® P, X = M x f2 and F=lfU {9}. Apply the same reasoning as in the proof of the 
second assertion in the last lemma, however, now with changing 9' , to get 

(5.2) ta(J8«' P )fj(i I u 1 k)6XxG : f#, r (a;) • d(T h , r (ar, w),T(a;, w)) > e}) = 0. 

Let H C G be those h for which 

Um ^({(a:,w) &K : p% ir (x)d{T h , r (x, u), T(x, u)) > e}) > 0. 



Because we know that (5.2) holds, we must have linv^oo 9' r {H) = 0. Hence, there are countably many 
g £ G such that 

lim 6 (Ux,u) G K : d(T 9tr (x,w),T(x, u)) > el) = 0, 

where we used that Pg r /*• 1 for (A* ® P) a.e. (x,uj), according to the Lemma above. This shows that 
the Proposition holds for those g. 

Pick one such g £ G. Then the first part of the previous lemma implies 

Q g B r -» Q°° vaguely on M X M X Cl. 

implies that for ai 

by Corollary 



This in turn impl ies that for any h £ G {Th)*Q g B r ~ > {i~h)*Q co Q°° by invariance of Q°° . Moreover, 
we have (Th)*Q g B r — QhgB r - This means, that for any h £ G we have 



3.12 



QhgB T Q°° vaguely on M x M x Q. 
Applying once more the first part of the previous Lemma proves the Proposition. □ 
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Corollary 5.10. There is a measurable map <J : M(M) x M(M) -> M(M x M) s.t. <p := *(A",^) 
denotes the unique optimal semicoupling between A" and p u . In particular the optimal semicoupling is a 
factor. 

Proof. We showed that the optimal semicoupling Q°° can be constructed as the unique limit point of a 
sequence of deterministic functions of A* and p* . Hence, the map u) i— > q u is measurable with respect to the 
sigma algebra generated by A* and p,'. Thus, there is a measurable map iff such that q' = iff(\',p,'). □ 

5.2.1. Semicouplings of A* and a point process. If p' is known to be a point process the above convergence 
result can be significantly improved. Just as in Theorem 4.8 and Corollary 4.9 of |HS10j we get 

Theorem 5.11. For any g € G and every bounded Borel set A C M 

lim (A* ®P) ({(x,uj) e A x Q : T g r (x, w) ^ T(x, u)}) = 0. 

i — >oo 

Corollary 5.12. There exists a subsequence (r{)i such that 

T gri (x,uj) — > T(x,uS) as I — >• oo 

/or almost every x € M , cu £ and every g € G. Indeed, the sequence (T g ^ ri )i is finally stationary. That 
is, there exists a random variable l g : M X fi — > N sucft i/iai almost surely 

Tg jTl (x,u)) = T(x,u>) for all l>l g (x,u). 

6. The other semicouplings 

In the previous sections we studied semicouplings between two equivariant random measures A* and p* 
with intensities 1 and /3 < 1 respectively. In this section we want to remark on the case that p* has 
intensity fj > 1. Then, g* is a semicoupling between A* and fi* iff for all to G f7 

(tti),^ = A w and (7r 2 ),^ < 

This will complete the picture of semicouplings with one marginal being absolutely continuous. In the 
terminology of section [2.2| we should better talk about semicouplings between \x* and A*. However, we 
prefer to keep A* as first marginal as it better suits our intuition of transporting a continuous quantity 
somewhere. We will only prove the key technical lemma, existence and uniqueness of optimal semicou- 
plings on bounded sets. From that result one can deduce following the reasoning of the previous sections 
the respective results on existence and uniqueness for optimal semicouplings. We will not give the proofs 
because they are completely the same or become easier as we do not have to worry about densities. 

Lemma 6.1. Let p G L l (M,m) be a nonnegative density. Let p be an arbitrary measure on M with 
fi(M) > (p ■ m){M). Then, there is a unique semicoupling q between (p ■ m) and p minimizing Cost(-). 
Moreover, q — (id,T)#(p ■ m) for some measurable cyclically monotone map T. 

Proof. The existence of one Cost minimizing semicoupling q goes along the same lines as for example in 
Lemma |3.2| Let q± be one such minimizer. As q% is minimizing it has to be an optimal coupling between 
its marginals. Therefore, it is induced by a map, that is q\ = (id, Ti)* (p ■ m). Let (?2 = (id, T 2 )„(/9 • m) be 
another minimizer. Then, q^ = h (qi + 92) is minimizing as well. Hence, (73 = (id,T^) if (p ■ m). However, 



just as in the proof of Lemma 3.2 this implies T± = T 2 (pm) almost everywhere and therefore q\ — qi. □ 



7. Cost estimate for Compound Poisson processes 

In this section we state some cost estimates for the transport between the Lebesgue and a 7— compound 
Poisson process. We consider M = R d ,A* = L the Lebesgue measure and p* a 7— compound Poisson 
process of intensity one with iid weights (Xi)i^,X\ ~ 7 In [HS10 a general technique was developed 
which allows to deduce upper estimates on the transportation cost by upper moment estimates of the 
random variable p'(A). In short, having good bounds on moments and inverse moments of p* allows to 
deduce transportation cost estimates. 
We have the following estimates on LP— cost 

Proposition 7.1. i) Let p > 1 be such that E[Xf] — 00 and i?(r) > r^ p ~ 1 ' <1 , then = 00. 

ii) Assume d < 2 and E[Xi] < 00. Then for any concave i) : [l,oo) — > K dominating 1) 

°° *(r) , 

iii) Assume d > 3 and p < po — 1 with 2 < p = sup{q : E[Xf] < 00} < 00. Then for $(r) = r p we 
have Coo < 00. 
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iv) Assume d > 3 and K[Xf] < oo for all p > 0. Then for t9(r) = r p we have < oo for any p > 0. 

Proof, ad i) Cqo can easily be bounded from below by the cost of transporting mass X± optimally into a 
single point. This transportation cost behaves for fixed X\ = X\(yl) like 

rcX]> d 

i?(r)r d - 1 dr > X\. 



\Z - a\> ap}<2- exp(-a(2 + p - 2y/l + p)) < 2 ■ exp ( -a 



o 

Taking expectation wrt X\ yields the desired result. 

The other claims are straightforward adaptations of the techniques from [HSlOj . We omit the details. □ 

Marko and Timar |MTllj constructed an allocation of Lebesgue measure to a Poisson point process in 
dimensions d > 3 with optimal tail behavior. In our language this means that there is a constant k 
such that the optimal mean transportation between a Lebesgue measure and a Poisson point process 
with cost function c(x,y) = exp(/c|a; — y\ d ) are finite in dimensions d > 3. Up to the constant n this is 
optimal (e.g. see USUI i. Their construction is based on an algorithm by Ajtai, Komlos and Tusnady 
AKT84] and uses two key properties of the Poisson point process, independence on disjoint sets and 
exponential concentration around the mean in big cubes. In the case of a 7— compound Poisson process 
the independence is inherited from the Poisson process. If we take 7 to be the exponential distribution 
one can show 

Lemma 7.2. Let Z = X)j=i -^i with N a Poisson random variable with mean a and (Xi)^^ a sequence 
of iid exponentially distributed random variables with mean 1 independent of N . For any < p < 1 it 
holds that 

4 8 

Hence, by using the very same algorithm as in [MTllj one gets for 7 the exponential distribution with 
mean one 

Proposition 7.3. Let d > 3 and p* a "/—compound Poisson process. Then there is constant k > such 
that optimal mean transportation cost between C and p' for the cost function c(x,y) = exp(K • \x — y\ d ) 
is finite. 

8. Stability 

As an application of the previous results, especially the existence and uniqueness results, we want to 
study stability properties of the optimal coupling between two random measures. Moreover, we will show 
some metric properties of the mean transportation cost. 

Given sequences of random measures (A* )„ £ n, (/i") n gN and their optimal couplings (?*)neN we want to 
understand which kind of convergence A* — > A* , /1* — > p* implies the convergence q* n — » q* , where q' 
denotes the/an optimal coupling between A* and p° . If for all n A„, X,p n ,p are probability measures q n 
the optimal coupling between A„ and p n (all transportation cost involved bounded by some constant) 
and A„ — > X, p n — > p weakly, then, by the classical theory (see Theorem 5.20 in |Vil09j ) . also along a 
subsequence q n —¥ q weakly, where q is an optimal coupling between A and p. 

A naive approach to our problem would be to ask for A* A- A* and p' n A p'. However, in this case let A* 

and p* be two independent Poisson point process and set A* = p* n = A*. Then, we indeed have A* A- X* 

and p* n -4 p* . Yet, the optimal couplings (q*) n eN, which are just q^(dx,dy) = 5 x {dy)X b ^ l {dx), do not 
converge to any coupling between A* and p* in any reasonable sense. Moreover, the couplings (<z*)„eN do 
'converge' (they are all the same) to some coupling q* with marginals being A* and A* having the same 
distribution as A* and p* . 

The next best guess, instead of vague convergence in distribution is vague convergence on M x M x O. 
Together with some integrability condition this will be the answer if the cost of the couplings converge. 

For two random measure A*, p* with intensity one and c(x, y) = d p (x, y) with p € [1, 00) write 



,/x*)= inf <£(q')= inf E 
<?*en es (A',At*) (j'en es (A',/i') 



/ d p {x,y) q'(dx,dy) 

JMxBn 



IMxBo 

We want to establish a triangle inequality for W p and therefore restrict to L p cost functions. We could 
also extend this to more general cost functions by using Orlicz type norms as developed in [Stullbj . 
However, to keep notations simple we stick to this case. 
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In this chapter, we will assume that all pairs of random measures considered will be equivariant and 
modeled on the same probability space (0,21, P). As usual P is assumed to be stationary. Moreover, we 
will always assume without explicitly mentioning it that the mean transportation cost is finite. 



Recall the disintegration Theorem |2.6| This will allow us to use the gluing lemma. 
Proposition 8.1. Let /i*,A*,£* be three equivariant random measures of unit intensity. 

i) Wp(A , ,Ai*) = & A w = ^" P-a-s- 

ii) W p (A*,/i*) = W P (/AA # ). 

iii) W p (A%/0<W p (A',r)+W p (r,M'). 

Proof, i) W P (A*, /i*) = iff there is a coupling of A* and /i* which is entirely concentrated on the diagonal 
almost surely, that is iff A w = ^ P— almost surely. 

ii) Let q* be an optimal coupling between A* and \x* . By definition 



For g,h € G put 



W?(AV) = 

f(g,h)=E 



E 



d p {x,y)q'(dx,dy) 



MxBq 



d p (x,y)q*{dx,dy) . 

I gB xhB 

By equivariance and stationarity, we have f{g, h) — f(kg, kh) for all k € G. Hence, we can apply the 
mass transport principle. 



E 



h£G 



d p (x,y)q'{dx,dy) 



gB x M 



J2f(9,h)=E\[ 
sec IJa 



d p {x,y)q'(dx,dy) 



MxhBo 



This proves the symmetry. 

iii) The random measures are random variables on some Polish space. Therefore, we can use the gluing 
Lemma (cf. |Dud02j or |Vil09j . chapter 1) to construct an equivariant random measure q' onMxMxM 
such that 

(7Ti,2)*g* € n opt (\',n') and (7r 2 ,3)*?* G Ilop^/i*,^*), 
where nop^A*,^*) denotes the set of all optimal couplings between A* and /i*. q' is equivariant as 
the optimal couplings are equivariant and q* is glued together along the common marginal of these two 
couplings. 

To be more precise let q° G n opt (A*, fi') and q' G n opt (/i*, ^*). Then, consider lMxgB xQ<l' and 
lgB xMxnq* to produce with the usual gluing Lemma a measure q* on M x M x M x O with the 
desired marginals on M x g£> x M x fl. As all these sets are disjoint we can add up the different q* 
yielding q' = J2 g eG 1g a measure with the desired properties. 
For g,h e G put 

e(g,h)=E / d p (x,z)q'(dx,dy,dz) . 

J MxgB xhB 

By equivariance of q* , we have e(kg 7 kh) = e(g,h) for all k € G. By the mass transport principle this 
implies 



E / d p (x 1 z)q°(dx,dy,dz) =E 

J MxB xM 

Then we can conclude, using the Minkowski inequality 



MxMxBq 



d p (x, z)q'{dx, dy, dz) 



/ P (A*,n 

< E 
= E 

< E 



d p (x, z)q*(dx, dy, dz) 



M x M x B 



d p (x, z)q*(dx, dy, dz) 



M x Br, x M 



1/p 



1/p 



1 1/P 



d p (x, y)q m (dx, dy, dz) 

I M x Bo x M 

= w p (\',^) + w p (fi',C)- 

In the last step we used the symmetry shown in part ii) 



E 



d p (y, z)q'(dx, dy, dz) 



M x Sn x M 



1/p 



a 
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Remark 8.2. Note that the first two properties also hold for general cost functions and general semi- 
couplings. The assumption of equal intensity is not needed for these statements. 

Let V p = {equivariant random measures \x* : W p (m, /i*) < oo}. 

Proposition 8.3. Let (/i*)„ e N,M* € V p be random measures of intensity one. Let g* denote the optimal 
coupling between m and /i* and q* the optimal coupling between m and fj,' . Consider the following 
statements. 

i) W p (/i* , (J,*) -> as n -> oo. 

ii) /i*P — > /x'P vaguely and W p (/i* , m) — > Wp(/i*, m) as n — > oo. 

iii) g*P — > g*P vaguely and W p (/i* , m) — > W p (^*, m) as n — > oo. 

iv) g*P — » g*P vaguely and 



lim lim sup E 

-R->°o n->oo 



(C(B )r)xB 



<P{x,y)q* n {dx,dy) 



where (Bq)r denotes the R— neighbourhood of Bq. 
Then i) implies ii). iii) and iv) are equivalent and either of them implies i). 

Proof, i) ii) : For any / s C C (M x f2) we have to show that limn^oo E[/j„(/) — /i(/)] = 0. To this end, fix 
/ G C C (M x fi) such that supp(/) C K x f2 for some compact set if. / is uniformly continuous. Let rj > 
be arbitrary and set e = n/(2m(K)). Then, there is 6 such that d(x, y) < S implies d(f(x, ui), f(y, uj)) < e. 
Put A = {{x,y) : d(x,y) > 6} H M x K and denote by k* an optimal coupling between p? n and //*. By 
assumption, there is JV G N such that for all n > TV we have W£(/x*,/i*) < M/ife^TO ■ Then, we can 
estimate for n > N 



< 



E 



(f(x,u) - f(y,w))K%(dx,dy) 



MxM 



< e-m(K) + 



E 



(f(x,u) - f(y,w))K%(dx,dy) 



< 
< 



2||/|UE[<(A)] 



1 



2||/|| 0o ^(M;,/* , )-m(J0 



< V , »7 



2 2 

The second assertion in ii) is a direct consequence of the triangle inequality: 

W p «,m) <W P 0CO+W p (M',m) 



and 



Taking limits yields the claim. 



,0i',m) < W p (Ai* , a*') + W p (aC m). 



ui) iv) : By the existence and uniqueness result we know that q^(dx,dy) = 5T^( x ){dy)m(dx) and 



5.9 



q u (dx,dy) — 8 T „ ^(dy)m(dx) . In particular, we have [j% (dx)¥(du) — (T^)^m(dx)V(doj). By Lemma 
we know that the vague convergence of g*P — > g*P implies that T n —> T locally inm®P measure. This 
in turn implies the convergence of / o (id, T n ) — > / o (id, T) in m ® P measure for any continuous and 
compactly supported function / : M x M —> M. Then, it follows as in the proof of Lemma 5.9 that 



E J f(x,T n (x))m(dx) — > E / f{x,T(x))m(dx). 

Let Ck(x, y) be a continuous compactly supported function such that for any (x, y) € (i?o)fc-i X B Q we have 
d p (x, y) = Ck(x, y), for any x £ C(B )k we have Ck(x, y) = and ct(x, y) < d p (x, y) for all (x, y) G M x M. 
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Then, wc have 



lim sup E 

n— > oo 

< lim sup ( E 

7 

= E 



C((B ) R )xB 

d p {x,y)q' l (dx,dy) 



MxB 

d p (x,y)q'(dx,dy) 



M x S 



< E 



d p {x,y)q' n {dx,dy) 

E 1 
d p (x,y)q*(dx,dy) 



CR(x,y)q^(dx,dy) 



MxBq 

cn(x,y)q'(dx,dy) 



MxB a 



/C((Bo)ji-ixS 

Taking the limit of R — > oo proves the implication Hi) iv). The other direction is similar. 

iv) => i) : We will show that W p (/i*,/i*) — > by constructing a not optimal coupling between ^t* and 
/i* whose transportation cost converges to zero. Let T„,T be the transportation maps from the previous 
steps. Put Q n (dx, dy) := (T n ,T)*m. This is an equivariant coupling of /i* and fi* because the maps T„, T 
are equivariant in the sense that (see also Example 2.11) 



T e ^(x) = gT"(g~ 1 x). 



The transportation cost are given by 

C(On)=E 



dP(x,y)Q n (dx,dy) =E / d" (T n (x), T(x))m(dx) 

We want to divide the integral into four parts. Put A R — {x : d{T{x) 1 x) > R} and similarly A R = {x : 
d(T n (x),x) > R}. The four parts will be the integrals over B Q n E a A% n CM fl with a,b £ {0, 1} and 
C°A = A. We estimate the different integrals separately. 



E 



B nCA R nCA R 



d p (T n (x),T(x))m(dx) 



0. 



by a similar argument as in the previous step due to the convergence of T n —¥ T locally in m®P measure 
and the boundedness of the integrand. 



E 



B nA R nA R 



d p (T n (x),T{x))m(dx) 



< 2 P E 



B n nA R 



d p (x,T(x))m(dx) 



+ 2 P E 


/ d p (x,T n (x))m(dx) 




JB nA R 



If d(x,y) < R,d(x,z) > R and d(y,z) < d(x,z) + R + a for some constant a(= diam(_Bo)), there is a 
constant Ci, e.g. C\ = 2+diam(f? ), such that d(y, z) < Cid(x,z) (because d(x, z)+R+a < (2+a)d(x, z)). 
This allows to estimate with (x = x,T(x) — z,T n (x) = y) 



Similarly 



E 



E 



B Q nCA R nA R 



d p {T n (x),T(x))m(dx) 



I B nA R nCA R 
This finally gives 



d p (T n {x),T(x))m(dx) 



<C P E 



< C p E 



B nA R 



d p (x,T(x))m(dx) 



B nA R 



d p (x, T n (x))m(dx) 



lim sup E 



d p (x,y)Q n (dx,dy) 



B x M 



< lim lim sup 2 P E 



Cf E 



B nA R 



B nA R 
d p (x,T(x))m(dx) 



d p (x,T(x))m(dx) 
C p E 



2 P E 



B nA R 



' B Q nA R 
d p (x, T n (x))m(dx) 



d p (x, T n (x))m(dx) 



= 0, 
by assumption. 



□ 
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Remark 8.4. For an equivalence of all statements we would need that ii) implies iii). In the classical 
theory this is precisely the stability result (Theorem 5.20 in |Vil09| ) . This result is proven by using 
the characterization of optimal transports by cyclical monotone supports. However, as mentioned in the 
discussion on local optimality (see Remark |4.4[ ) a cyclical monotone support is not sufficient for optimality 
in our case. 

We do not have real stability in general but we get at least close to it. 

Proposition 8.5. Let (A* ) n <=n and (/i*) n eN be two sequences of equivariant random measures. Let g* be 
the unique optimal coupling between A* and /i*. Assume that A*P — > A*P vaguely, /i*P — > /i*P vaguely 
and sup„ £(g*) < c < oo. Then, there is an equivariant coupling g* of A* and fi' and a subsequence 
(g* ,)fceN such that g* P — > g*P vaguely, the support of q* is cyclically monotone and 

£(<?*) < liminf £(g*). 



In particular, if 



hm C(0 = . 



g* is the/an optimal coupling between A* and /x* and g*P — > g*P vaguely. 

The proof is basically the same as for Proposition |2.20| Hence, we omit the details. 

Remark 8.6. The last proposition also holds if we consider semicouplings instead of couplings (see 



Proposition |2.20[ ) 



Example 8.7 (Wiener mosaic). Let fx' be a Poisson point process of intensity one on K 3 . Let each atom 
of no evolve according to independent Brownian motions for some time t. The resulting discrete random 
measure is again a Poisson point process, denoted by \i° t (e.g. see page 404 of |Doo53| ). Consider the 
transport problem between the Lebesgue measure C and /i* with cost function c(x,y) = \x — y\ 2 . Let q* 
be the unique optimal coupling between C and /i*. Then, (£((?*) = W%(C, /tx*) = W2 (£,£**) for any seM 
as n' s and /z* are both Poisson point processes of intensity one. Moreover, we clearly have /i*P — > /x*P 
vaguely as s — > t and therefore q*P — > g*P vaguely. By Lemma [5. 9| this implies the convergence of the 
transport maps T s — ¥ T t locally in C ® P measure. In particular, we get a continuously moving mosaic. 

Example 8.8 (Voronoi tessellation). Let /i* be a simple point process of unit intensity. Put Hg=j3- fi*. 
We want to consider semicouplings between the Lebesgue measure C and /i* for fi > 1. By the results 
of section [6j there is a unique optimal semicoupling q' between C and /x*. Moreover, q^ = (id,Tp)*£ 
for some measurable map Tp : R d x SI — > Mr. It is clear that g^, induces the Voronoi tessellation with 
respect to the support of ^ no matter which cost function ■& we consider. We want to show that g* — > g* 
vaguely as (3 — > 7 for large 7. For this it is sufficient to show that C : (3 H> £(g*) is continuous. For 
7 > P > 1 note that g^ is also a semicoupling of £ and Hence, C(/3) is monotonously decreasing in 
j3. Moreover, from the previous Lemma we know that 

<£(?;)< liminf 

Therefore, C(/?) is a monotonously decreasing lower semicontinuous function. This implies that it is right 
continuous. With a bit of work it is also possible to show that C(-) is left continuous in f3 < 00. For 
f3 = 00 one can show directly that g* — > g^, vaguely. 
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